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FOREWORD 


Lockheed Missiles & Space Company, Inc., Huntsville Research & 
Engineering Center, conducted the effort reported herein under Contract 
NAS8-29824 for the National Aeronautics and Space Administration, 
George C. Marshall Space Flight. Center. 

The work was performed during the period May 1973 through August 
1974. The NASA contracting officer's representatives for this contract 
are Mr. James W. Bilbro and Mr. Harold B. Jeffreys, alternate. 
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Section 1 

INTRODUCTION AND SUMMARY 

In May 1973, Lockheed-Huntsville received a contract from NASA- 
MSFC to: (1) develop a servo-controlled scanner for NASA's Laser Doppler 
Velocimeter (LDV); (2) to perform certain field tests of LDV components; 
and (3) to develop a mathematical model of the scanning LDV interacting 
with aircraft- generated vortices. In December 1973, the contract was ex- 
panded to cover additional testing, certain hardware modifications and com- 
plete scanner manual development. Again in June 1974, the contract was 
expanded to include mechanical service testing at the MSFC test site and to 
perform planning for operation at a remote site. 

Lockheed developed three range and elevation scanners — one set for 
each of three MSFC LDV units — and installed these subsystems at MSFC in 
late 1973 and early 1974. The elevation scanners provide a capability to 
manually point the LDV telescope at operator chosen angles from 3.2 to 89.6 
degrees within 0.2 degrees or to automatically scan the units between operator 
chosen limits at operator chosen rates of 0.1 to 0.5 Hz. Likewise, the range 
scanners provide a capability to manually adjust the focal point of the system 
from a range of 32 meters to a range of 896 meters under operator control 
or to scan between operator chosen limits and. at rates from 0.1 to 6.9 Hz. 

The range and elevation scanner controls are designed to allow simultaneous 
range and elevation scanning so as to provide finger scan patterns, arc scan 
patterns and vertical line scan patterns. The development and testing of these 
units are discussed in Section 2 of this volume and also in Volume II of this final 
report. 

Numerous field and laboratory tests have been performed under this 
contract to evaluate vendor-provided components, to optimize marginal sub- 
systems and to overcome unforseen hardware difficulties. Lockheed personnel 
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have tested numerous optical components and subassemblies of the LDV 
both in the field and in the laboratory. Lasers, interferometers and specific 
optical elements have been examined both independent of the system and 
while assembled into the system. Calculations and experiments have been 
performed to determine system efficiency from experimental and theoretical 
standpoints. Experiments were conducted along with analytical efforts to 
develop a practical means of reducing backscattered radiation from the 
secondary mirror of the telescope which previously saturated the detector 
and reduced the effectiveness of the Mach Zehnder interferometer. The 
translator, which earlier had been judged very inefficient, was shown to have 
promise of efficient operation with modifications. Spatial resolution of the 
system was measured both in the laboratory and in the field. A longer-than- 
expected spatial resolution was attributed to laser beam truncation by the 
telescope primary mirror. This truncation was also found to reduce system 
signal-to-noise ratio from hard targets at long ranges due to the increased 
size of the focal volume. However, system efficiency appeared quite adequate 
to detect wake vortices during numerous aircraft flybys at MSEC. Laboratory 
and field tests are discussed in Section 3 of this volume. 

Lockheed developed a fluid dynamic model of the wake vortex and 
modeled the interaction of the laser system with the vortex. The fluid dynamic 
model of the vortex included such effects as the centrifugal force upon the par- 
ticulate matter suspended in the vortex. The LDV model included the system 
spatial resolution and time of passage of particles through the focal volume. 
Details of this modeling effort are presented in Section 4. 
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Section 2 

SCANNER development AND TESTING 

Authors: E. W. Coffey 

B. B. Edwards 

Lockheed-Huntsville received a contract to fabricate three combination 
range -elevation scanners for the NASA-MSFC Laser Doppler System in May 
1973. The units were delivered the following November, December and 
February, The units consisted of independent range and elevation servo 
systems with a digital control panel for each. The elevation scanners were 
designed to be pointed, using a manual mode, from 3.2 to 89 deg accurate to 
within 0.2 deg and scan in an automatic mode within operator chosen limits 
and at operator chosen rates from 0.1 to 0.5 Hz. The range scanners were 
likewise designed to be pointed manually from distances from 32 to 599 meters 
and scanned automatically between operator chosen limits and at rates from 
0.1 to 6.9 Hz. Various operator designated inputs would program the units 
to perform finger scans, arc scans or coordinated line scans. 

The details of scanner design and operator instructions are presented 
in Volume II, "Scanner Operations Manual," of this final report. However, 
a description of results of the acceptance tests plus additional tests per- 
formed on the scanners is presented in the following paragraphs. Section 
2.1 presents the results of the acceptance test on unit 001. Section 2.2 pre- 
sents the results of additional tests which were used to optimize the range 
calibration technique. 

Tests referred to in this section as E-1, E-2, etc., refer to tests run 
to evaluate the "elevation" scanning capability of the unit. Tests referred to 
as R-1, R-2, etc., refer to tests run to evaluate the range scanning capability 
of the unit. 


3 



LMSC-HREC TR D390159-I 


2. 1 RESULTS OF ACCEPTANCE TEST ON LASER SCANNER 001 

A description of the tests agreed upon as being suitable is given in the 
Appendix. This report presents the results obtained when these tests were 
carried out on unit 001, along with some comment as considered appropriate. 


Test E" 1 

In this test the output of the D/A converter, which operates on the ele- 
vation command angle 12- bit word, was recorded on a strip chart with a known 
calibration. The front panel thumbwheel switch, which controls the upper 
limit of the auto- scan mode was advanced successively in 1-deg increments 
from 10 to 59 deg. The purpose of this test was to verify that the command 
word did, in fact, reverse direction when it reached an angle corresponding 
to the dialed- in upper limit, over the full range of variation of this variable. 
The other conditions of this test were: 

a. Mode: auto- scan, normal 

b. Lower limit: 0 deg 

c. Scan frequency: 0.1 Hz 

The results are shown in Fig. Z-L Note that the turnaround angles 
correspond exactly to those of Table A- 1 of the Appendix. Since the computer, 
which senses when the upper limit has been reached, uses only the five most 
significant bits of the word it is evident that on the ^*up*’ cycle the command 
word corresponding to equality on the comparator for the first time will have 
all its lower order bits (7 through 1) equal to zero. All possible dialed-in 
upper limit angles having identical structure in bits IZ through 8 will thus 
produce the same effect on the command word, despite the fact that they 
represent different angles. This truncation effect creates the "plateaus" 
in the input/output table of Fig. Z-1. 

The system has since been modified to reduce truncation effects. 


4 
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Test E-2 

This test was identical to E-1 with the roles of upper and lower limit 
interchanged. The results are shown in Fig. 2-2. As ejjpected, complete 
conformity to the input/output characteristic of Table A- 2 of the Appendix 
was obtained. (We observe that the command word on the "down" angle will, 
when the e<juality condition on bits 12 through 8 is met for the first time, 
have all the lower order bits (7 through 1) in the "high" condition, creating 
the "staircase" effect of Table A-2. The actual scan will always be "inside" 
the dialed-in limit values.) 

Test E- 3 

In this test the scanner was made to execute the normal auto- scan 
function in elevation between six conditions of upper and lower limits as 
shown in Table A-3, with five (the complete set) selected scan frequencies 
in each distinct upper/lower limit condition. The purpose was to verify that 
the logic, which is supposed to vary the counter clock rate to achieve the 
same dialed-in scan frequency for any condition of limit boundaries (i. e., for 
any scan excursion) was, in fact, accomplishing this. The recorded command- 
word cycle behavior is shown in Figs. 2-3-1 through 2-3-6. The scan fre- 
quency was computed from the known chart speed, with results as shown in 
Table 2-1. Note that all the frequencies have a slight bias toward high values. 
This arises from the truncation effect mentioned in Tests E-1 and E-2. The 
control logic uses the entire 12-bit upper and lower limit words (and the 4-bit 
BCD scan frequency control word) to compute and generate the appropriate 
counter clock rate, but the scan turn-around occurs automatically between 
the limits set by the 5-bit comparators. Since this excursion is always 
"inside" the true one (i.e., smaller), the rate is slightly high. It was agreed 
that this effect was acceptably small. (The bias-in rate could, if desired, be 
split (^) to have the correct average value by a minor hard-wire change in the 
fixed rate multiplier, which is now set at lZ/64 ratio.) 
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Fig *2-3-1 - Scan Frequency Test (Elevation Scanner Between 58 and 41 deg) 
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1.5 deg 



Fig. 2-3-2 - Scan Frequency Test (Elevation Scanner Between 42 and 25 deg) 
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Table 2-1 

MEASURED SCAN RATES OF TEST E3 


Angle 

Designation 

Excursion Limits 
(deg) 

Scan Rates (Hz) 

0.1 

0.2 

0.3 

0.4 

0.5 

A 

58-41 

0.1065 

0.216 

0.32 

0.430 

0.534 

B 

42-25 

0.107 

0.211 

0.32 

0.434 

0.527 

C 

26-8 

0.113 

0.224 

0.337 

0.445 

0.565 

D 

58-25 

0.104 

0.206 

0.315 

0.415 

0.515 

E 

42-8 

0.107 ' 

0.213 

0.32 

0.427 

0.54 

F 

58-8 

0.107 

0.208 

0.306 

0.419 

0.524 
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Test E-4 

The purpose of this test was to ascertain that the anti-crossover logic 
was operating correctly. This logic is constructed so as to substitute, 
for the dialed- in lower scan limit, a new limit equal to the (actual) upper 
limit minus 9.6 deg, whenever the difference between the upper limit and the 
dialed- in lower limit is less than 9.6 deg (and particularly whenever the dialed- 
in lower limit is larger than the upper limit, i. e. , crosses over). The test 
consisted of setting the upper limit at three different values and, for each 
value, moving the lower limit in successive steps toward and beyond the 
(fixed) upper limit while recording the actual scan excursion. Figures Z-4-1 
and 2-4-2 show the results, i.e., the lower limit approaching the 0^^ - 9.6 deg 
value* and staying there, thus demonstrating the correct operation. 

Test E- 5 


The purpose of this test was to demonstrate correct operation of the 
manual mode slewing function and the single increment capability. This was 
done by: 

a. Operating both the fast and slow slewing through the full range 
and recording the command voltage to show the slew rate and 
automatic top and bottom limiting, and 

b. Depressing the "single increment',' button a large number of 
times in succession in both the up and down directions and 
recording the successive step increments in the command 
word on a calibrated scale to show the step size. 


The results are shown in Figs. 2-5-1 and 2-5-2. The slew rates were 
within 10% of their design values (34 deg/sec and 3.4 deg/sec, and automatically 
limited at 89 deg and 3.2 deg as required. 


In certain instances the 5-bit truncation effect is in evidence on - 9.6 deg. 
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40-29 deg 







Fig. 2-5-2 - Incremental Step Size Check 
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The single increment operation occasionally yielded steps greater than 
the design value of 0.2 deg, indicating that the one-shot gate controlling this 
operation is slightly longer than it should be, but was considered acceptable. 

Test E-6 

The purpose of this test was to demonstrate the proper operation of the 
sector scan mode. In this mode the scan limits for the automatic mode are 
derived from a pair of 8-bit counters. One of these counters controls the 
spread between the upper and lower limits, moving both symmetrically toward 
or away from a center value, while the second moves the center value up or 
down across the elevation domain- 

While the controls exhibited all the correct functions, it was clear that 
the slew rates used in the design were too fast, and this test was deferred 
pending a decision on proper slew rates. The rates can easily be reduced 
by a factor of 10 by minor modifications to the card containing the clock gen- 
eration circuitry. 

Tests R-1 Through R-6 

These tests were identical to their elevation counterparts, mutatis 
mutandis. Since the circuitry is identical, acceptable performance by one 
confirms acceptable performance by the other. 

The following series of tests was made to ascertain the static and 
dynamic accuracy of the servo systems. 


Slew rates have been decreased. 

The range scale factor in meters is ten times that of the elevation variable 
in degrees, and the scan frequency range spans from 0.1 to 6.9 Hz, requiring 
two BCD control words. 


17 
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Test E-7-1 

In this test, made on the elevation servo (manual mode), a series of 
consecutive 5 - deg increments was set with the manual slewing control, 
and the actual elevation mirror angle was read from the vernier dial mounted 
on the mirror rotation axis. The percent accuracy is shown in Table 2-2. 
Some comments are in order; 

• A zero offset was not properly nulled out 
so that percent error is abnormally large 
at the small angles. 

• The percent error exceeds 0.2 at all points 
tested. 

• A significant amount of hysteresis is exhibited. 


Most of the error observed may be attributed to insufficient "stiffness” in the 
servo to overcome friction and small inertial moment unbalance. The system 
at present has a static error factor of 0.123 deg/in-oz of friction or mass un- 
balance. Although the static error coefficient could be reduced by the use of 
a higher gain in the error channel (Gj of Fig. 2-6), there are two reasons why 
this is not possible: 


1. The use of higher gain Gj, necessitates increasing gain G 3 
in the damping (tachometer feedback) channel in order to 
avoid overshoot transients at the turn-around points of the 
sawtooth driving waveform. If both these gains are made 
high, the other transient term in the response, which is 
due to the real root of the third order polynomial of the 
transfer function, reaches excessive duration and magnitude. 

2 . The dynamic range at the final power amplifier (which sub- 
tracts the tachometer voltage from the error voltage) is 
exceeded, and the resultant response is very much under- 
damped, since proportionality cannot be maintained when 
saturation is reached. 


The LED readouts were used to determine the command angle actually arrived 
at with the slewing control. 
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Table 2-2 

STATIC POSITIONING ACCURACY IN ELEVATION ANGLE 


Command 

Angle 

(deg) 

Specified 

Lower 

Limit 

Specified 

Upper 

Limit 

Actual 

Percent 

Error 

Command 

Angle 

(deg) 

Specified 

Lower 

Limit 

Specifi ed 
Upper 
Limit 

Actual 

Percent 

Error 

5 

4.975 

5.025 

4.95 

-1 

65 

64.675 

65.325 

65,4 

+0.615 

10 

9.95 


10.05 

+0.5 

70 

69.65 

70.35 

70.35 

+0.5 

15 

14.925 


15.0 

0 

75 

74.625 

75.375 

75.0 

0 

20 

. . 

19.9 


20.1 

+0.5 






25 

24.875 

25.125 

25.2 

+0.8 






30 

29.85 1 

30.15 

30.15 

+0.5 






35 

i 

34.825 

35.175 

35.25 

+0.714 






40 

39.8 


40.5 

+ 1.25 






45 

44.775 

45.225 

45.3 

1 

+0.666 






50 

49.75 

50.25 

50.4 

+0.80 






55 

54.725 

55.725 

55.2 

+0.364 






60 

59.70 

60.3 

60.15 

+0.25 
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For these reasons any significant reduction in the static error can only be 
obtained by a more exact balance of the large mirror mass distribution 
throughout the region of elevation angle being scanned. It should be noted 
that these errors, provided they are not too large, are not of major impor- 
tance, since the actual pointing angle of the mirror is read by the feedback 
potentiometer (which is accurate to 0.1% of full scale) and small discrepancies 
between the command position and the actual position are tolerable (the com- 
mand angle is not part of the recorded data). 

Test E-7-2 

During this test, simultaneous strip chart, recordings were made of the 
command angle and the actual elevation mirror angle with the system in the 
auto-scan mode. Several combinations of scan excursion and scan rate were 
examined, the purpose being to determine the phase lag at high rates, the 
reduction in scan excursion, and the '^rounding off” of the response at the 
turn-around points. The results are shown in Figs. 2-7-1 through 2-7-20. 

It is clear from the response illustrated in these figures that the servo 
bandwidth is undersirably low. Unfortunately, this cannot be remedied by 
increasing the loop gain since the constraint is the upper bound on the available 
amplifier voltage driving the torque motor. It can readily be shown that 
the maximum frequency to which the servo can respond,. If the maximum 
acceleration is applied at all times is 

f = 1/4 (a (2.1) 

max ' max' s' ' ' 


here 


and 


a 

max 

e 

s 


2 

maximum aval la ble acceleration (rad/sec ) 
total scan angle (rad) 


a 

max 



E 

max 


/R J 


(2-2) 
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Fig. 2-7-5 - 0 to 59 Degree Excursion at 0.1 






















Fig. 2-7-11 - 45 to 59 Degree Excursion 
at 0. 1 Hz 



Fig. 2-7-12 - 45 to 59 Degree Excursion 
at 0.2 Hz 
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Fig. 2-7-15 - 45 to 59 Degree Excursion 
at 0.5 Hz 



Fig. 2-7-16 - 10 to 25 Degree Excursion 
at 0. 1 Hz 
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Command 



Fig. 2-7-19 - 10 to 25 Degree Excursion 
at 0.4 Hz 


Fig. 2-7-20 - 10 to 25 Degree Excursion 
at 0.5 Hz 
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where 


= motor torque constant (Ln-oz/amp) 

E a maximum available amplifier voltage (volts) 

R = torque motor resistance (ohms) 

2 

J = moment of inertia of load plus motor (in-oz/sec ) 


In the actual design the fixed parameters are; 


K,p = 60 in-oz/amp R a 9 ohms 

Emax ” volts J - 38 in.-oz/sec^ 


which gives 


Vax = rad/sec' 


For an angular excursion at 0.5 rad (28.6 deg) this yields 


f 

max 


0.7 Hz. 


Since this value corresponds to a condition at 90-deg phase lag and continuous 
acceleration, the actual limiting frequency will depend on the distortion criteria 
imposed — it certainly must be smaller than that given above. A substantial in- 
crease in bandwidth would be possible if the torque motor had a much lower 
resistance. 

Test R-8 

This test involving range scanning against an external wheel target at a 
known range, was performed later and is discussed in Section 2.2 of this report. 


Arising from the saturation, i.e., full torque is always applied and only the 
direction is changed to create the periodic reversals of direction. 
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Test R-9 

Figure 2-9* shows the display x, y coordinates for successive 50 meter 
range increments; the dotted lines correspond to inward displacement and the 
solid lines correspond to outward displacement. The following was noted: 

(1) A hysteresis error, increasing linearly to about 8 meters at 500 meters, 
and (2) Jitter on the x, y voltage due to the truncation effect of the digital (13- 
bit optical encoder) feedback. It was decided to defer the x, y accuracy tests 
until more work had been done in adjusting offsets and gains in the range 
channel. It was also discovered that the feedback gain was set to a value equal 
to one-half of its correct value; this was later corrected. 

Test R-10 — Range Servo Dynamics 

The dynamic response for a variety of input range commands was 
examined by simultaneous comparison of the input range command and the 
reconstructed range voltage. The results are shown in Figs. 2-10-1 through 
2-10-12. The main effect at high scan rates is to reduce slightly the over- 
all excursion in range, and to introduce a very small phase lag. (Note: The 
actual mirror position is obtained from the feedback voltage, so that the out- 
put X, y does not exhibit any lag.) 

Test R- 1 1 — Line Scan Mode 

This mode was tested by setting in several values of horizontal range, 
programming the scanner through an angle at 60 deg scan and observing the 
X, y on an x, y recorder. Ideally, the x, y plots should be vertical lines 
through the x-axis values corresponding to the horizontal range. We note 
some hysteresis effect in the results. This test should be repeated since 
later realignment of the x, y channel obtained much better accuracy. 

To maintain consistency of test and figure numbers. Fig. 8 is intentionally 
omitted. 

**In the figures, the input command is the lower trace and the reconstructed 
range voltage is the upper trace. 
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Fig, 2-10-3 - 30 to 500 Meter Fig. 2-10-4 - 30 to 500 Meter Excursion 

Excursion at 6. 9 Hz at 3 Hz 
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-10-6 - 30 to 200 Meter Excursion at 1 Hz 
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Fig. 2- 10-7 - 30 to 200 Meter Fig, 2-10-8 - 30 to 200 Meter Excursion at 6.9 Hz 

Excursion at 3.0 Hz 
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5 , 2-10-10 - 300 to 500 Meter Excursion 
at 1 .0 Hz 
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Fig. 2-10-11 - 300 to 500 Meter Excursion FI 
at 3.0 Hz 




Lg. 2-10-12 - 300 to 500 Meter Excursion 
at 6.9 Hz 
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2.Z RESULTS OF ADDITIONAL SCANNER TESTS 
2.2.1 Introduction 


The data from the previous flyby tests indicated the following: 


« The range Coverage of both systems was considerably smaller 
than the programmed value. 

• The system in van 2 gave elevation angles consistently lower 
than the programmed values, at least on the data readouts. 
(The actual angle, however, as judged visually, appeared 
correct. ) 

Investigation of these discrepancies revealed the following: 


9 Note from Fig. 2-11 that the scale factor relating secondary 
mirror position to input range command is determined by 
the ratio, G^/G 2 - adjusting this to the correct value, a 
value of 

Gg = IO/Ztt volts/radian 

was used, based upon full scale voltage (lOV) corresponding 
to full rotation (0 = Ztt radians). Since only the lower 12 bits 
from the optical encoder is actually used for feedback, the 
true gain is 

G_ = IO/tt radian/volt 

D 

i.e. , Gc is really twice as large as was thought. The system 
was aligned for secondary mirror position offset by placing a 
target at a known distance 


Rj =61 meters (200 feet), 

programming the input command to this value, and displacing 
the secondary mirror platform to give optimum focus on this 
target. With the incorrect value of G j/g^ used, this resulted 
in the relation 


R 


actual 


R 


command 


1/2 + 


R 


command 


2R 


1 
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We observe that the actual range agrees with the command 
only when the command is equal to Rj, but that for other 
values the scan coverage will be reduced. This was corrected 
by readjusting Gj/G 2 to the correct value. 

The discrepancy in angle was found to exist only on the readout data for 
van 2, and was due to the method used to correct the sine/cosine potentiometer 
for the left-handed configuration required in van 2. This resulted in a dis- 
tortion of the sine function. The condition was corrected by using the same 
potentiometer connections as in van 1, and inverting the sine voltage with a 
unity gain operational amplifier (the previous configuration used this amplifier 
in the non -inverting mode). 

After these corrections were made, both systems were realigned and 
calibrated with respect to the x, y outputs. This procedure involves adjusting 
a series of gains and offsets in the chain of devices generating x and y. This 
data output chain is shown in Fig. 2-12. (The y output is the same, with the ex- 
ception that the cosine function replaces the sine function. ) After alignment, 
the values 


6 - = tan"^ (V /v ) 
el ' y X 


„ , 2. 2| l/2 

R = K • V + V " 
six y 


were compared with the actual elevation angle (as read on the vernier dial on 
the elevation mirror shaft) and with the actual external focus range as computed 
from the relation between the optical encoder output voltage and the secondary 
mirror position. (This comparison relies on the validity of the lens equation 
relating secondary mirror position to external focus, and on the precision with 
which the optical encoder output indicates the true mirror position. ) 


Results of the calibration are shown in Figs. 2-13 and 2-14. We note that 
angle errors are less than 0.5 degree and range errors are a small fraction of 
the external range resolution length. 


In van 2, the opposite (from van 1) rotational sense on the potentiometer shaft 
occurs to indicate increasing (i. e. , upward) elevation angle. This is because 
van 2 views the other end of the scan plane baseline, but employs the same 
directional layout of optical axis. 
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Command 

Angle 

(deg) 


Command 

Range 

(tn) 

Actual 

Range 

(m) 

tan’^ (y/x) 
(deg) 

, 2. 

(x + y ) 

(m) 



50 

50.386 

5.67 

49.44 

5.3 

5.3 

100 

99.85 

5.31 

97.65 



200 

199.0 

4.92 

194.35 



350 

352.6 

5.03 

341.51 



50 

50.41 

15.52 

48.47 

15.0 

15.2 

100 

99.76 

15.20 

97.65 



200 

201.40 

15.27 

197.71 



350 

352.09 

15.14 

344.87 



50 

50.38 

30.66 

49.91 

30.0 

30.3 

100 

100.32 

30.43 

99.66 



200 

198.6 

30.39 

197.43 



350 

352.09 

30.25 

345.7 



50 

50.66 

45.36 

50.66 

45.0 

45.1 

100 

99.86 

45.17 

100.36 



200 

200.09 

45.05 

200.20 



350 

340.7 

45.16 

346.30 


Fig. 2-13- Output Data Static Accuracy for Van 1 System 
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Fig. 2 -14a - Errors in R = (x + y ) ' After Calibration for Van 2 
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In addition to the above, tests were made on the shrinkage of the range 
scan excursion as a function of scan rate (an inevitable consequence of a finite 
bandwidth servo with max torque limitations, i. e. , any physically realizable 
servo, being required to follow commands involving extremely large accelera- 
tions) and the hysteresis error between command and output range. These are 
shown in Figs. Z-I5 and 2-16. Neither of these effects is serious, since the data 
readout is correct and the programmed scan excursion can always be made 
to yield any desired actual scan coverage; in addition, the range excursion 
is quite linear, hence there is no difficulty with variation of integration time 
in the processor. 

2.2.2 Tests on Hard Target (Wheels) 

After calibration the system was scanned against wheel targets at 
several known ranges. To obtain an indication of the behavior of the external 
focus, the signal back-scattered from the wheel was used as the input to a 
spectrum analyzer set in the non-scanning mode, with its IF bandwidth cen- 
tered on the Doppler frequency, and with a video bandwidth sufficient to pass 
the expected modulation frequencies without distortion. Under these circum- 
stances we would expect the signal (which is proportional to back scattered 
power) to indicate the variation in intensity incident upon the wheel, and to 
exhibit a peak at a range equal to the known target range. To determine if 
this was actually the case, the spectrum analyzer output and the range voltage 
were simultaneously recorded on a Brush recorder. Figures 2-17 and 2-18 show 
the results obtained at two different ranges, the computed range (obtained by 
projecting the time of occurrence of peak signal on to the calibrated range 
sweep) agrees fairly closely with the known range value. 

There are two difficulties in using such a procedure: 

1. The peak signal does not necessarily occur in time coincidence 
with the incidence of maximum intensity, since the effective 
backs catter is the sum of a large number of randomly phased 
vectors reflected from individual small grains of particulate 
matter (sandpaper) on the wheel surface. 
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Fig. 2-15 - Range 
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Input 

Command 

Outward 
T raverse 

^■act 

Inward 
T raverse 

\ct 

50 

49.66 

50.39 

100 

99.38 

100.435 

150 

147.25 

150.07 

200 

194.78 

199.96 

250 

241.67 

249.56 

300 

289.24 

298.51 

350 

333.44 

344.41 


Notes: 1. Lack of symmetry with respect to input indicates 
gain or offset error or both. 

2. Actual R determined from encoder output. 


Fig. 2-16 - Static Hysteresis — Van 2 Range Servo 
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Fig* 2-17 - Signal Returned from Short Range Target 
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-3 dB 
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Signal from Belt Sander . 

Video Bandwidth = 100 Hz 
Scan Freq. = 0.1 Hz 
Range = 229 m 
Run 1, Van 1, s/ 11/1 A 


Fig. 2-18 - Signal Returned from Long Range Target 
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2. When the offset on the secondary mirror position is removed 
by using a target at short range, any errors in finding this 
mirror location are greatly magnified when the mirror is 
moved to a position corresponding to long range. 

While it is true that the axial resolution of the system also is magnified 
in the same proportion as is any setting error, this does not remove the error, 
but merely makes it to be the same percentage of the axial resolution. It would, 
then, be best to perform the offset alignment with a calibration target at long 
range, using a very long averaging time to determine the setting corresponding 
to peak value. Any setting errors resulting would be then quite small at short 
range. 

Since the strip chart data give the variation of signal level with focus 
range, it should then be possible to determine the range interval between 3dB 
down points, i. e. , the axial resolution. This was done for the ranges involved, 
with results as shown in Fig. 2-19. The results are comparable to these obtained 
previously on the same type of optical system using direct measurements of 
the intensity variation. In Fig. 2-19 'mode' refers to heterodyning configuration. 

2.2.3 Air Flow Targets 

To determine the ability of the system to locate an air flow of velocity 
higher than the ambient air, an airplane propeller was set up at a range of 152 
meters (500 feet) and the system scanned against this in several modes. Fig- 
\ires 2-ZO, 2-21, and 2-22 show scans corresponding to: 

• Fig. 2-20 — Range only, with angle such as to pierce flow stream 

• Fig. 2-21 — Angle only, with range fixed at 152 meters 

• Fig. 2-22 — Range and angle, both dimensions enclosing target. 

The expected variations are all present, i. e. , 

• The peaking of the signal on the range scan, with minima 
at the extremities of the scan (Fig. 2-20). 

' The recording method was as described in the preceding section, except 
that now the IF filter was moved up to the neighborhood of 1.5 MHz, so 
as to not respond to the (continually present) wind background signal. 


53 



LMSC-HREC TR D390159-I 


Wheel, 250 ft 


Scan Rate (Hz) 

0.1 

1 

2 


6 

Mode 

Van 

Resolution (ft) 

Front Cavity 

1 

15.7 

21.86 

24.49 

26.24 

36.7 


2 

X 

24.20 

23.70 

25.80 

X 

Secondary Lr. O. 

1 

13.3 

34.98 

34.98 

21.50 

34.9 


2 

X 

31.16 

35.50 

38.80 

X 

Separate L. 0. 

1 

16.9 

mgm 

X 

28.90 

36.2 


2 

28.9 

Bi 

33.90 

32.30 

X 


Belt, 750 ft 


Scan Rate (Hz) 

0.1 

1 

2 

4 

6 

Mode 

Van 

Resolution (ft) 

Front Cavity 

1 

162.7 

1 

150.6 

X 

141.5 

X 

Secondary L. O. 

1 

136.8 

137,8 

X 

144.0 

X 

Separate L. Oa 

1 

134.9 

144.6 

X 

142.7 

X 


Fig. 2-19 - Signal Returned from Long Range Target 
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Fig* 2 -20a - Signal Returned from Airstream - Slow Scan (Range Only) 
















Signal from Blower 

Video Bandwidth = 100 Hz (fixed bandpass spectrum analyzer) 
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Fig. 2 -20b - Signal Returned from Airstream - Fast Scan (Range Only) 












Fig. 2-Zl - Signal Returned from Airstream (Angle - Only Scan) 
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Fig» 2-22 - Signal Returned from Airstream (Range and Angle Scan) 
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• The (more or less) continuous signal available over a finite 
angular region, with noise only outside this region, as ex- 
hibited by angle -only scan (Fig* 2-21). 

• The group, finite in number, of periodic signal maxima, 
with "dead space" between such clusters, exhibited by the 
range and angle scan. 

Two computations were made from these data: (1) the average range to 

the peak signal; and (2) the physical extent of the "signal-present" region in 
the angular scan -only mode. The range measurement average gave a location 
of 170 meters, i.e., about 18 meters longer than the correct value (the in- 
dividual readings were always within + 3 meters of this average value). The 
second measurement gave a value of 16 feet for the diameter of the airstream 
containing velocities significantly higher than the wind. This compares quite 
well with the known diameter of the propeller span, which was 14 feet (some 
spreading out of the airstream is to be expected). 

The range discrepancy was later checked against a wheel target, and 
was found to be still present. The system was then adjusted for zero offset 
against a target at 700 feet, and the location of the peak determined for a wheel 
at 500 feet. The location thus determined was within +4 meters of the true 
location. This result emphasizes the importance of a very careful initial 
offset alignment, preferably done at long range . It also emphasizes that 
the major error in locating a flow field will not be due to inaccuracies in the 
system readouts, but in the displacement of the peak signal due to random 
signal fluctuations arising from the relative locations and backscatter magni- 
tudes at the aerosol particles comprising the target. In order to locate accu- 
rately, the intensity pattern (illumination) must move across the particle 
stream at a speed such that the same assemblage of particles is present dur- 
ing the entire time taken to change the illumination by a significant amount. 

If a new array of particles enters the beam while the focus is moving, it will 
be impossible to distinguish signal level changes due to the changing illumi- 
nation from signal level changes arising from the new array of scatterers. 

The latter can easily create 10 dB or more signal change. 


'^This is difficult to do, since movement of personnel in the van can result in 
a slight tilt of the floor, changing the location of the incident spot on the wheel 
and the corresponding Doppler shift. In addition, fluctuations at signal level 
are always present, obscuring the location of the maximum. 
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A series of similar tests scanning against the same blower airstream, 
now moved to 85 meters range, was made several days later. Some of the 
results, recorded in the manner just described, are shown in Figs, 2-23, 2-24 
and 2-2 5. While the general nature of the signals is similar to those previously 
obtained, we note that only at the higher scan speeds was it possible to obtain 
a consistent estimate of the target location. One would expect that at short 
range the slower scan speeds would be acceptable, if we consider the following: 


Let 


Tp = particle population renewal time 

T = 

p V sin 0 

P 


T = scan modulation time 
s 

_ A (R/d)^ 

S R 

Tp/Xg = (R/vp sin e ) x d/R 


where 

X 

R 

d 

V sin 0 
P 

R 


wave length 

range to airstream 

effective optical aperture diameter 

cross-beam velocity of particles 

scan velocity 


We note that other things being equal, T /t varies as 1 /r. It is also 

/ 4 P ® 

true that at short range, due to the 1/R behavior of the illuminated volume, 
the ratio of mean to standard deviation of the wind background signal is very 
much diminished (the volume of the higher speed particles in the airstream 
is constant) so that the occasional large particles at wind velocity may cause 
excitation of the higher frequency filters (since their spectrum is spread out 
because of short signal duration) and obscure the scan modulation of the air- 
stream. Additional tests should be made to determine the quantitative nature 
of these effects. 
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Fig. 2-23 - Range and Angle Scanning {64 to 192 m at 4 Hz, 3.2° to 16° at 0.1 Hz) 
Typical Signal Segment - Propellor at R =85 m, Van 2, 11 June 1974 
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Fig, 2-24 - Angle-Only Scan (R = 85m, Target Range; 

0 = 3.2° to 16° at 0.1 Hz, Van 2, 11 June 1974) 
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a. 4 Hz Range Scan 

Fig. 2-25 - Range Scan, Propeller 
R = 64 to 192 m, Van 2, 
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In addition to the data presented in Figs. Z-20, 2-21, and 2-22 which were 
()l)tained with the use of a spectrum analyzer, the outputs from the processor 
were used to obtain estimates of target range and spatial resolution, under 
the condition of range-only scan with the airstream target at 152 meters 
range. To obtain these estimates, the signal strength as a function of range 
was plotted for several scans. The measure of signal strength used was the 
mean over 10 successive 4 msec integration intervals of the amplitude of the 
filter which displayed the peak indication (the wind peak being ignored). The 
data are shown in Fig. 2-26. 

The propeller wake was at a range of 152 meters but appears in the 
data as being at a range of 170 meters due to a scanner range misadjustment. 
Estimates of the separation between the 3 dB gain points of the curve yields 
a value (corrected for the range misadjustment error) ~ 25 meters. If the 
spatial extent of the flow field (~ 5 meters) is subtracted an estimate of the 
spatial resolution of ~ 20 meters at a range of 150 meters is obtained. This 
is in excellent agreement with the values obtained using hard targets obtained 
during the Astrionics testing. 
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Section 3 

LASER DOPPLER SYSTEM OPTICS TESTING AND OPTIMIZATION 

Authors: Dr. C-C. Huang 
L. K, Morrison 
J. L. Jetton 


Lockheed-Huntsville personnel have performed numerous laboratory 
and field tests on the optical components of the laser Doppler system under 
the subject contract. Laboratory optical tests were performed and docu- 
mented by C-C. Huang with L. K. Morrison and J. L. Jetton assisting. Various 
field tests were performed by each of the above individuals. The primary 
objectives of these tests were to: (1) evaluate Lockheed, Raytheon and other 

vendor-provided components of the laser Doppler systems; (2) optimize per- 
formance of marginal components; and (3) develop new techniques to surmount 
unforeseen difficulties (such as secondary backscatter)* 

Section 3. 1 includes discussions of the backscatter problem and several 
approaches which were tried to develop a solution to this problem. Initial 
Lockheed-Huntsville system calibration tests are discussed in Section 3.2. 
Spatial resolution and system signal-to-noise ratio tests conducted m the 
laboratory are presented in Section 3.3. Section 3.4 presents signal-to-noise 
ratio and spatial resolution measurements made in the field and compares 
these results with theoretical results. Section 3.5 presents tests results for 
the interferometer tests and polarization tests of the lasers used in the system. 
Section 3.6 provides a discussion of ideas for developing bias circuitry. The 
Isomet acousto-optical translator is discussed in Section 3. 7, 
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rcUirns, before reduction, make operation in tlie Mach Zehnder interferometer 
ctjiifigu ration impractical with a low power coherent detector. Without reducing 
tlie backscattered radiation from the secondary, this backscattered portion of 
the local oscillator power is larger than the desired local oscillator power. 

Methods suggested to reduce this undesirable local oscillator level 
have included use of: 

• Hole in 45-degree mirror 

• Spatial filter 

• Wire in front of secondary, and 

• Location of X/4 plate and lens between 
secondary and primary 

The analyses of these techniques are discussed in Section 3. 1. 1, while actual 
tests results of the first three techniques are presented in Section 3.1.2. 

3.1.1 Analysis of Backscatter Reduction Techniques 

• Hole in 45-Degree Mirror Approach 

A method proposed to reduce the backscatter from the secondary mirror 

of the LDV is shown in the sketch on top of the next page. The central part 

of the laser beam will be a "dark" shadow after reflection by the 45 -degree 

mirror with a 1.41 mm hole in the middle. Unfortunately the shadow will 

2 

close up at a Rayleigh distance (d /X) which is 10 cm away from the 45-degree 
mirror. The backscatter from the secondary mirror possesses a central 
"bright" shadow beyond a 10 cm distance. 

The diffraction pattern of the 45-degree mirror with a 1.41 mm hole 
can be produced analogously with aid of a CO^ laser beam incident upon a 
6 mm circular aperture with an opaque ball (1 mm in size). We can there- 
fore investigate the diffraction pattern of that ball as shown on the sketch at. 
the bottom of the next page. 
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These arrangements can be scaled longitudinally and laterally by a 
factor of 16,8 for easier diagonistics by utilizing an He-Ne laser rather than 
a CO^ laser. 

The following arrangement for the 0.6328 jjim He-Ne laser should lead 
to conclusions that can be scaled for CO^ lasers. 




Instead of a 0.36mm circular opening with an opaque ball (0,06 mm in 
diameter), a circular opening (0.4mm) with a metal wire (gauge No. 40 or 
0.003 in. or 0,076 mm in diameter) was used in the experiment to measure 
the diffraction pattern of the wire without loss of generality. 

When the screen was located just behind the aperture, the dark shadow 
of the wire was observed as expected. When the screen was placed 1 cm 
away from the aperture the central dark shadow could not be observed, A 
central "bright'* shadow of diameter ~0.2 mm was observed on the screen 
6 cm away from the aperture and the ''bright" shadow extended toward 
infinity. 
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We conclude that employing a 45-degree mirror with a hole in the middle 
is probably not the solution to eliminate the backscatter from the secondary 
mirror. 

• Spatial Filter Approach 

Another method proposed to reduce the backscatter from the secondary 
mirror of the LDV is shown schematically by the following sketch. 



The spatial filter will transmit collimated radiation through the pin 
hole and block the backscatter from the secondary mirror. In this way 
the backscatter can be reduced. However, as the actual reduction in back- 
scatter power is much less than 13 dB under most circumstances as ex- 
plained in the following analysis. 

Since the diffraction limitation can only focus a 6 mm CO 2 laser beam 
to a diameter = 0. 02Z6cm by using a 10 cm focal length lens, a 
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pin-hole size > 0.0226 cm is required for the spatial filter. The back- 
scatter from the secondary mirror can be regarded as a point source 60 cm 
away from the 10cm focal length lens and focused to a point 12 cm behind the 
lens. 

In the following figure we have: 

0 = X/ttCO , 20 =» 0 . 6/10 rad 

ooo 

2(0 = 4\/ZirQ 0,0226cm 

Ooo 

where X = 10.6/im. 



Filter 
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The spot size of the backscatter from the secondary mirror at the spatial 
filter [)osiiion can be expressed as 


2co(z) = 2(0 


T 1/2 


1 + 


where z = 2.0 cm in this case and gives 2co(z) = 3.84 x (20i^) - 0.1025 cm. 

The attenuation of the secondary mirror backscatter is much less than 
lOlog (203 {z)/2oo^^)^ = 10 log 21.0 = 13.23 dB because of the Gaussian distri- 
bution in laser intensity. 


The maximum attenuation of the backscatter from the secondary mirror 
has been calculated and plotted as shown in the following sketch versus the 
focal length of the lens used to focus the signal from infinity through a best 
matched spatial filter. 



Secondary Mirror 0*ck»cattor Attenuation Lrns 1 - oral Lcnmh 


*Kogelnit, H., and T. Lii> "Laser Beana and Resonators," Proc.IEEE, Vol54, 
No, 10, October 1966, pp. 1315-1316. 
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• Shifted X/A Plate and Lens Approach 

In the following paragraphs, a method (untested) is discussed of pre- 
venting undesired laser power being seen by the detector, the origin being 
the reflection of laser power from the telescope secondary mirror through 
the central hole of the telescope primary mirror. 

A diagram of this isolation technique is depicted below; 


Target 
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Most of the S“wave from the CO^ laser is reflected by the Ge coupler. 
The transmitted CO^ laser is approximately 100% p-wave. The laser is 
focused through the central hole of a X/4 plate by a lens and then reflected 
by the secondary mirror. The laser beam through the X/ 4 plate is then cir- 
cularly polarized and transmitted to targets via the primary mirror. The 
backscattered laser radiation from target is collected by the primary mirror 
and passed through the X/4 plate. After this stage the radiation from target 
is changed from circularly polarized radiation into a s-wave. After the col- 
lection by the secondary mirror through the central hole of the X/4 plate, this 
s-wave is reflected by the Ge coupler, mirror, beamsplitter and into the 
detector. 

The reflection loss from the central portion of the secondary mirror 
will retain the same polarization as the incident laser beam and will transmit 
through the Ge coupler with a polarization that will not be sensed by the 
detector. We can have nearly 100% isolation between the reflection loss from 
the secondary mirror and the aerosol signal only under the assumption that 
the outgoing laser beam through the Ge coupler is 100% linearly polarized. 

The configuration of the optics for this system is shown on the following 

page. 


The laser power loss due to this configuration can be estimated by as- 
suming that the intensity distribution of the outgoing laser beam from the 
telescope secondary mirror is gaussian over the X/4 plate. The laser power 
transmitted through the X/4 plate and focused upon the target is 


u 


P = P, 


o VF / 


max 


e du 


u 


min 
u 


= P 


O hf 


max Z ™ 

e du - p 


I 


^min 2 


e"^ du 
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Lens: CdTe Meniscus lens with diameter 1 cm and focal length +4" AR coated on both 

surfaces at i0.6/[jtm; Clear aperture > 7 mm, surface figure better than X/20 at 
10.6/ljtm over 80% area, 

QdS — Plate: Diameter = 1" jfO.Ol", Thickness <2 mm, diameter of the central hole = 3 mm 

' +0.1 mm. Edges chamfered < . 5 mm, AR coated on both surfaces at 10. 6/im, 

surface figure better than X/ZO at 10.6/J.m over 90% area. 


Secondary Mirror: Diameter = l/2" parabolic onvex Au coated mirror, f - -1.00 cm, surface 

figure better than X/20 at lO.bjUm over 80% area. 
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= 0.9545 

■^o 



o 


0.70 p. 


where p^ is the total laser power before entry into the telescope. So we 
have about 30% loss in laser power for this geometrical configuration. 

The A/4 plate in this proposed system will change slightly the axial 
focal depth because of its finite thicjcness (d), finite index of refraction (n), 
and the f/Z configuration of optics. The image size (Ap) of a point source 
seen through a slab of dielectric material with thickness (d) and index of re- 
fraction (n) can be expressed as 


Ap = (1 


COS 1 V J 

) d - 

cos r 




cos i . d 
cos r ' n 


which is a section of a line in the direction normal to the slab, where i and 
r are incident and refracted angles respectively for extreme rays and sini/ 
sin r = n. 


For CdS {n = 2,29) in f/2 system the above equation can be written as 


Ap 


(1 


cos i . ^ 
cos r ^ n 


(1 - 0.975)^ 


0.011 d 


Ap< 0.022mm (= 0.000866") for d<2mm. 
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The additional stretch of the axial focal volume (Aq) beyond the diffraction 
limit due to spherical abberration effects introduced by the CdS \/4 plate can be 
written as: 

r 

where 

p 24" 
f ^ 24" 

Ap <0,000866" 

and therefore, 

Aq <1.8 X lO"^ forq»f. 


Aq is estimated to be less than 20 feet for a range of 1000 feet. This should 
be an endurable tolerance. Aq is much smaller at shorter range, (See sketch 
on the following page. 
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Since the laser beam passes through the A/4 plate divergently for the 
outgoing beam and convergently for incoming beam, a polarization loss is in- 
troduced. The polarization loss is the laser power contributed to the original 
direction of polarization after a round trip through the X/ 4 plate. This part 
of the power will not be collected by the detector. 


The time independent part of the linearly polarized laser beam can be 
written as 


E 


1 
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After a round trip through the \/4 plate the extreme rays can be written 
as 


E 

ex 


i . 7T / ^d* 
"J 9 t A' 




\ 


(¥) 


E 


y 


where Ad is the optical path difference between the extreme rays and 
paraxial rays; in our case 


‘ = 0 : 97 - 1 =1-03 -1 = 0.03 for £/2 system. 


therefore 


• 0. 01 Stt 


\ 


ex 


J • O.OIStt 


j 


The electric field which passes through a polarizer with its axis along the 
original polarization direction is 
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The laser power passing through this linear polarizer is 


(E ) • E 
P P 




2 sin 

V2 

2 sin 0.01 5;^ 


= sin 0.01 Stt 
= E^ sin 2. 7° 

= 0.047 E^ 


This is a rough estimate but we can see that, due to the X/4 plate, the 
polarization loss of the backscattered signal should be less than 5%. 


3.1.2 Secondary Mirror Backscatter Reduction Experiments 


• 45-Degree Mirror with Hole for Backscatter Reduction 

The experimental set-up is shown on the sketch on the following page. 
A piece of asbestos paper was placed very close to the 45-degree mirror 
as depicted in the sketch. The diffraction pattern of the laser beam which 
was produced by the 45-degree mirror appeared on the asbestos paper as 
a bright doughnut shape. The central dark region filled in and turned 
bright when the asbestos paper was moved closer to the secondary mirror. 
Although the 45-degree mirror with the hole was used, the laser beam in- 
cident upon the secondary mirror appeared with a Gaussian-like intensity 
distribution as expected. However, the transmitted laser beam intensity 
was reduced by a factor of 6.6 {8.2 dB) because of the hole in the middle 
of the 45-degree mirror. 


A smaller attenuation factor should apply to the signal receiving 
leg because the backscattered power reflected by the 45-degree mirror with 
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the hole has an approximately uniform illumination of power upon the 45 -degree 
mirror instead of a Guassian beam attenuation in the transmission leg. 

The total signal strength is linearly proportional to the product of the 
signal power and the local oscillator power, but each was reduced by approxi- 
mately 8,2 dB because of the transmitted power loss through the hole of the 
45-degree mirror. A 16.4 dB drop in signal level was expected in addition 
to the secondary backscatter reduction by this arrangement. 

The experimental results are shown in the following two pictures on 
the next; page. The total reduction in signal level was 16 dB which was es- 
sentially the same reduction factor as due to the laser power drop without 
the backscatter reduction effect included (i.e., the system is probably shot 
noise limited). 
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UmODtJCIBILrrY OF THE 
ORIGINAL PAGE IS POOR 



a. 45-Degree Mirror With Hole, Signal Level: -62 dBm 



b. 45-Degree Mirror Without Hole, Signal Level: -46 dBm 

Sandpaper Wheel Signal and Secondary Mirror Backscatter 
Reduction by a 45-Degree Mirror With Hole 
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A clear conclusion can be drawn, based on experimental and theoretical 
evidence, that the 45-degree mirror with the hole cannot reduce the backscatter 
power efficiently. 


• Spatial Filter and Backscatter Reduction 


The experimental set-up is shown below. 



The spatial filter consists of two identical lenses with a focal length of 2" and 
a 105|Hm pinhole. These three pieces were aligned carefully to allow maximum 
secondary mirror backscatter power through the spatial filter. The 105jUm 
pinhole could be translated along the optical axis as well as in two directions 
both perpendicular to the optical axis. The distance between the secondary 
mirror and the pinhole was 200 cm. 

While the signal measurement was taken, the pinhole had been adjusted 
along the optical axis to peak up the signal reading. The result is shown in 
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picture (a) on the following page. When the pinhole was removed, the signal 
was increased by 2.5dB in picture (b). Another 1 dB increase in signal level 
appeared when the spatial filter was totally removed as shown in picture (c). 
These results show that the backscatter reduction using this particular spatial 
filter set-up was Z,5dB at most. 

• Liamp-Blacked Wire and Backscatter Reduction 
The experimental set-up is shown below: 
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is ?Oqj 


a. Wheel Signal With the Spatial 
Filter Set-Up and Pinhole 

Signal: -35.5 dBm 



b. Wheel Signal With the Spatial 
Filter Set-Up but Without 
Pinhole 

Signal: -33 dBm 



c. Wheel Signal Without Spatial 
Filter and Pinhole 

Signal: -32 dBm 


Sandpaper Wheel Signal and Secondary Mirror Backscatter 
Reduction by a Spatial Filter 
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With the separate local oscillator (L.O.) leg blocked, the wheel signal in the 
backscatter mode was recorded as shown in picture (a) on the the following 
page. A 0,4mm lamp*blacked wire was arranged to move across the front 
surface of the secondary mirror carefully to a position such that the mini- 
mum signal was recorded as shown in picture (b). The signal dropped by 
16 dB, mostly due to the reduction of secondary backscatter, A small per- 
centage (10%) reduction of the transmitted laser power onto the samdpaper 
wheel due to the wire was observed. 

When the separate L.O. leg was unblocked, the signal level reverted 
back to almost the same signal level as that of the backscatter mode without 
the wire (effect shown in picture (c)). The backscatter reduction by a 0.4 mm 
lamp-blacked wire was around 16 dB. Meanwhile the reduction of the signal 
was small. 

In the sketch shown below, the 0.4 mm lamp-blacked wire was 
vertically in front of secondary mirror and translated horizontally on the 
LDV table across the laser beam. The thermopile reading, which was 
directly proportional to the secondary mirror backscatter power (shown on 
the following page) is plotted as a function of wire position. When the wire 
approached the center of the secondary mirror, the edge diffraction pattern 
of the wire from the divergent backscatter power of the secondary mirror 
was clearly seen on the dc microvolt meter reading. 
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a. Backscatter Mode With Wire 
Out 

Signal; -32 dBm 


b. Backscatter Mode With Wire 
In 

Signal: -48 dBm 


c. Separate L.O. Mode With 
Wire In 

Signal; -30 dBm 


Sandpaper Wheel Signal and Secondary Mirror Backscatter 
Reduction by a 14 mm Lamp Blacked Wire 




88 




LMSC-HREC TR D390 159-1 


Before the wire reached the secondary mirror center, a very strong 
diffraction peak — almost twice as strong as the backscatter power without 
the wire — was measured by the thermopile. This strong diffraction peak 
from thd secondary mirror backscatter and the wire is potentially dangerous 
to a detector which cannot accept this high amount of power. As the wire 
was moved farther toward the secondary mirror center, a minimum amount 
of power - 16 dB smaller than the backscatter power with the wire removed 
- was measured by thermopile. Although the wire diameter was only 0.4 mm, 
the shadow of the wire will not be "filled in" due to the diffraction because of 
the large divergent backscatter power of the secondary mirror. Usually a 
small hump of the thermovoltage is recorded while the wire is placed right 
at the middle of the secondary mirror. This relates to the fact that the dif- 
fraction of a sharp-edged circular opaque disk placed in front o f a point source 
produces a bright spot along the axis .which connects the source and the 
center of the disk. Continued translation of the wire produces another mini- 
mum and finally the edge diffraction is seen again. As seen by the results 
depicted in the following curve, the small hump at the bottom part of the curve 
is missing due to the coarse incremental translation of the wire. The 16 dB 
attenuation of the secondary mirror backscatter power recorded by this meas- 
urement agrees very well with the coherent method. 

3.2 CALIBRATION OF THE LDV SYSTEM AT LOCKHEED 
3.2.1 Sulphur Wheel Target Calibration 

A 9. 5 watt CO^ laser beam was focused through the LDV telescope onto a 
rotating sulphur wheel located 100 feet from the telescope. The diffusely re- 
flected radiation from the wheel was collected by a 2-inch focal length BaF 2 
lens with a 1.9cm clear aperture, and focused onto a calibrated thermopile 
for the incoherent measurement. The distance between the laser spot on the 
wheel and the collecting lens was 100 cm. 

The theoretical value of the diffusely reflected power that the lens 
should receive can be expressed as shown in Eq. (1) on page 91. 
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MEASUREMENT OF THE RELATIVE BACKSCATTER POWER 
FROM SECONDARY MIRROR AS A FUNCTION OF WIRE POSITION 



Traveled Distance of a 0.4 mm Wire Across the Front of Secondary Mirror (mm) 
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P 


r 



P. cosO 
1 r 


U) 


for a Lambertian surface with reflectivity p, where 


A = clear aperture area of the collecting lens 

R = distance between the laser spot on the wheel 
and the collecting lens 

P. = incident laser power onto the sulphur wheel 

e = angle between the normal of the wheel and the 
^ direction of the reflection. 

For constant incident power and reflectivity (Eq. 1) the diffusely reflected 
power from the surface is independent of the incident angle but depends on 
the reflection angle as cos6^ . This has been proven experimentally and 
the resulting data tabulated in the table below. 


Laser 
Incident 
Angle (deg) 

B 

Experiment 

Theory 

Pj.,exp 
Pr, th 

hh 



BB 

Pj./ cos 6^ , mW 

IT ' 
(%) 

0 

33 


0.64 


0.810 

79 

45 

12 

IM 

0.60 

mm 

0.810 

74 

60 

27 

0.50 

0. 56 

0.722 

0.810 

69 

0 

33 

0. 54 

0.64 

0.679 

0.810 

79 

0 

48 

0.46 

0.69 

0.542 

0.810 

85 

0 

63 

0.31 

0.69 

0.368 

0.810 

85 

0 

73 

0.20 

0.69 

0.237 

0.810 

85 
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Within the jflO% error of the experiment (because the response of the 
thermopile was too sensitive to its position) the sulphur wheel surface can 
be referred to as a Lambertian surface with a reQectivity of 80% on average. 

The measured diffuse reflection of various targets as a function of the 
targets' orientation is shown by the curves on the following page. The exper- 
imental configuration is also shown. A stabilized CO^ laser was directed 
on the target when at an angle of 15 degrees to a line joining the target and 
the thermopile. The diffusely reflected power from the target that passed 
through a calibrated aperture was collected by a Ge lens and focused into 
the thermopile. The thermopile output voltage was a direct indication of the 
collected diffuse reflection power from the target. The target was positioned 
such that a line normal to its surface lay in the same plane with the incident 
laser beam and the thermopile. The target surface was continuously rotated 
to avoid local heating by the strong laser beam. The orientation of the target 
surface could be varied by rotating the target about an axis perpendicular to 
to the plane which includes the incident laser beam and the thermopile. 

The normalized reflectivity (normalized relative to the 0-degree inci- 
dence angle condition) is shown for three different targets; (1) sulphur wheel, 
(2) sandpaper wheel, and (3) belt sander. 

The experimental results show that the sulphur wheel provided a uni- 
form reflectivity over a wide range of incidence and reflection angles. This 
indicates that sulphur provides a good Lambertian surface. The absolute 
reflectivity of the sulphur wheel was measured to be 80% for the CO^ radiation 
as mentioned previously. 

The belt sander and the sandpaper wheel also can be regarded as a 
Lambertian surface for small incidence and reflection angles. The sandpaper 

with the smaller grain size provides a better match to a Lambertian surface 
then the coarser sandpaper. The absolute reflectivity of the sandpaper sur- 
face is approximately 8% 


92 



alized Reflectivity 



LMSC-HREC TR D390159- 





LMSC-HREC TR D390159-I 


3.2.2 Losses of the LDV System 
• Losses of Optics 


The sketch shown below gives the arrangement of the optics in the LDV 
set-up. 


T cli’ siopt' 

( ^ 



Sulphur 

Wheel 



Beamspl ilte r 


R =78.8%, 


T = 19.6% 


The loss to the outgoing laser beam was 42%, which did not include the 
first 45-degree beamsplitter transmission loss. {The beamsplitter was in- 
stalled to provide another leg of the local oscillator (L.O. ) power and can be 
replaced by a highly reflective 45-degree mirror.) 


94 




LMSC-HREC TR D390159-I 


The individual contributions to the total outgoing loss from each optical 
component were: 

Vertical Ge Brewster angle coupler loss: 14% 


S-wave component of the laser output: 0.24% 

45-degree mirror and \/4 plate loss: 10% 

Telescope loss: 20% 

Output mirror loss: 6,7% 


In addition to the spherical aberration and truncation effect, which will 
be discussed later, the incoming signal undergoes a total loss of 54% through 
the optics before reaching the detector. 

The individual contributions to the loss of the incoming signal due to 
the optics were: 


Output mirror loss: 6.7% 

Telescope loss: 20% 

45- degree mirror and X/4 plate loss: 10% 

S-wave coupling loss: 22% 

Mirror and beamsplitter loss: 5% 

Horizontal Ge Brewster angle coupler loss: 2% 

Collecting lens loss: 5% 


• Spherical Aberration and Truncation Loss 

Analysis of an f/2 system conducted by Lockheed-Huntsville shows that 
the spherical aberration contribution to the lateral size of the focal volume 
at a range of 100 feet and below and the spherical aberration contribution to 
the longitudinal size of the focal volume at a range of 250 feet and below 
should be considered seriously. A truncation effect exists at any range as 
long as the LDV telescope destroys the Gaussian nature of the transmitted 
waveform. 
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Both the lateral and the longitudinal elongation of the focal volume due 
to the spherical aberration and the truncation effect makes the signal focal 
volume {which was collected by the primary mirror) a prolonged oval shape. 
Since the virtual object of the secondary mirror was of a prolonged oval shape 
the signal beam collimated by the secondary mirror inherited a large divergent 
angle. The signal power at the detector position could not be fully recovered 
by using a 6 mm diameter aperture (or even 1.2 cm) in front of the collecting 
lens for any range. The measured results will be discussed later (Section 3.2.5h 

3.2.3 Optimization of PbSnTe Detector Bias Circuitry 

The bias circuit used in the LDV system calibration experiment has 
been modified for better impedance matching between the detector and the 
preamplifier. The modification of the circuit included adding a pulse trans- 
former, which has a bandwidth from 10 kHz to 5 MHz and a voltage gain of 
15.5dB, for the correction of the amplifier loading problem. Using the 
modified bias circuit, the highest signal-to-noise ratio recorded from the 
sulphur wheel at a 50-foot range was 104.5 dB under the operating condition 
of 10 mW L.O. power from the backscatter of the secondary mirror, zero 
bias voltage, 7.5-degree orientation of the wheel to the laser beam, and 10 W 
laser power incident on the sulphur wheel. 

The signal and the noise levels versus bias voltage are shown on the 
following page for various L.O. power levels for the Raytheon PbSnTe detector 
314, chip 1. It is obvious that the highest signal-to-noise ratio of that particular 
detector was obtained at zero bias voltage condition. 

3.2.4 Theoretical Signal -to -Noise Ratio of the Sulphur Wheel Doppler 
Signal and System Efficiency Tests 

The signal-to-noise ratio, S/n, for a reverse-biased photodiode is 
given by 


(S/N) 


power 


= T|Pj./hi f 


( 2 ) 
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O L>,0. power = 1.46 mW 
® L.O. power = 2.92 mW 


Sulphur wheel signal for a 100 ft range 
(213 cm between the secondary mirror 
and the 1.2 cm aperture) 


Q L.O. power = 10 mW 


Sulphur wheel signal for a 50 ft range 
(155 cm between the secondary mirror 
and the 1.2 cm aperture) 


Sulphur Wheel Signal and Noise Levels Under Various 
Operation Conditions 
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where 

T) = quantum efficiency of the detector ^ 15% 

-34 

h = 6.625 X 10 Joules-sec 
j/ = 3 X 10^^ Hz 
Af = 30 kHz 

Pj, ^ signal power incident upon the detector 

i 

therefore 

(s/n) =%* 2.5 X lO^^^P 

power r 

With the aid of Eq. (1), Eq. (3) can be rewritten as 


<S/%ower=“^ 5xl0‘^Tp ^P.cose^ (4) 

Typically, 

T = transmission of the LDV system to the return signal 
power ^46% (1-0,54) 

p = 80% 

Z 

A = area of the telescope primary mirror = ;r{15cm) 

R - range = 50 ft = 1500 cm 

p. = low 

X 

0 =45 degrees 

r ^ 

and we have that 


(S/N) 6.5x10 

'power 


10 
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or 

10 log (S/N) 108 dB 

® ' 'power 

For various calibration arrangements, Eq. (4) can be tabulated as 
follows; 


R 

(ft) 

6r 

(deg) 

P. 

(S/N)nower. dB 

Theory 

Experiment 

50 

45 

10 

108.0 

-- 

50 

7.5 

10 

109.5 

104.5 

100 

45 

10 

102.0 

95.5 

100 

7.5 

10 

103.5 

98.5 


The experimental value of the signal to noise ratio was about 5 dB below 
the theoretical value which can be fully explained in the following section. 

3.2.5 Verification of the Spherical Aberration and Truncation Effect 
• Incoherent Measurement 

The experimental set-up was essentially the same as shown in 
Section 3.2.2. The PbSnTe detector was replaced by a Sensors, Inc., ther- 
mopile to measure the return signal power. The distance between the 1.2 
cm aperture, which was located in front of the thermopile, and the secondary 
mirror was 213 cm. For a 50-foot range, the measured signal power was 
only about 0.017 mW when a piece of 0.4mm diameter wire was placed across 
the front of the secondary mirror to reduce the backscatter power. This 
enabled us to measure the slight difference in thermal voltage when the signal 
power was superposed upon the existing large backscatter power from the 
secondary mirror. 

As the thermopile was moved closer to the secondary mirror from 213 
cm to 141 cm, the measured signal power was 0.042 mW which was 2.5 times 
stronger. This proves the large divergence of the signal power due to the 
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spherical aberration and truncation effect. By referring to the table in Sec- 
tion 3.2.2 and the experimental arrangement (shown in Section 3.2.2) used 
to obtain the above result, we can estimate the signal power that should be 
received by the thermopile if there were no spherical aberration or trun- 
cation effect. Since the lens used in the experiment subtended the same 
solid angle as the telescope primary mirror for a 50-foot range, the total 
signal power received by the LDV system should be 0.46 mW. After a 20% 
two-way loss due to the wire and the 54% one-way transmission loss of the 
system, the signal power should be 0.17 mW which is four times higher than 
0.042 mW. 

• Coherent Measurement 

For the coherent measurement we maintained the distance from the 
secondary mirror to the aperture in front of the PbSnTe detector at 155 cm 
and measured the heterodyne signal from the sulphur wheel which was oriented 

7.5 degrees to the incident laser beam and placed 50 feet away from the tele- 
scope primary mirror. When the aperture was changed from the diameter 
of 1.2 cm to a diameter of 6mm, the signal- to- noise ratio dropped by 3.5dB. 
Only 1 dB drop of the signal-to-noise ratio was measured for the same pro- 
cess at a 100 foot range. This illustrated the larger divergence for the closer 
range due to the spherical aberration and the truncation effect. 

3.2.5 Conclusions 

The 104.5 dB measured signal-to-noise ratio off the sulphur wheel was 
only 5 dB below the theoretical value. This is due to spherical aberration 
and truncation effect. If there were no other additional losses of the signal 
power because of the range variation, we would expect less than 6 dB reduc- 
tion in signal-to-noise ratio as the range is changed from 50 feet to 100 feet 
solely because the spherical aberration decreases with increasing range. The 
observed reduction in signail-to- noise ratio was still 6 dB (refer to table on 
preceding page), which illustrated that the truncation loss is more dominant 
than the spherical aberration loss. 
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3.3 LDV SYSTEM SPATIAL RESOLUTION AND SIGNAL- TO- NOISE 
RATIO TESTS 

3.3.1 Test Conditions 


• A 13-inch folding mirror was placed inside a 100-meter 
tunnel to achieve the range of 100, 188, 3 50 and 558 feet. 

The 700-foot range was achieved when the folding mirror 
was placed at the outside end of the tunnel. 

• A sandpaper wheel was adapted as a hard target which 
provided a reflectivity of 7% when it was oriented 30 
degrees to the incident laser beam. 

• The incoherent power measurement at the nominal focal 
volume was made with various size pinholes attached in 
front of the power meter head. 

• Honeywell stabilized lasers were used during the test. 

• The interferometer was Mach-Zehnder Interferometer 1 
with modified mount on the detector which provided a three- 
dimensional translation and also with the modified mount on 
the recombining beamsplitter which provided rotations along 
two perpendicular axes and translation along the direction of 
the separate L.O. beam direction. 

® A 12-inch f/z telescope was used. 

• The air turbulence and the mechanical vibration of the LDV 
table degraded the accuracy of the data. 

3.3.2 Incoherent Measurements 
• Focused Spot Size vs Range 

The focused laser beam diameter at different ranges was determined 
with the help of a power meter and a pinhole attached in front. Both were 
moved across the focused beam at the focal plane. The total distance traveled 



power reading, was recorded as the diameter of the focused spot size. 
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Range 

(ft) 

Pinhole Size 
(mm) 

Focused Spot Size 
(mm) 

100 

0.5 

2.1 

188 

0.5 

4.5 

350 

1.0 

7.5 

558 

1.8 

12.7 

700 

1.8 

13.8 


These values are plotted in Fig. 3-1. 

# Focal Depth vs Range 

The focal depth at different ranges was determined by moving the power 
meter head, with a pinhole attached, along the optical axis at the nominal focal 
volume location. Then the distance span between two points was measured 
where the power meter reading dropped to half of its peak reading. 


Due to the lack of a longer optical bench, no attempt was made to measure 
the focal depth at ranges 558 and 700 feet. 


Range 

Pinhole Size 

Focal Depth 

(ft) 

(mm) 

(ft) 

100 

0.5 

2.9 

188 

1.0 

10.5 

350 

1.0 

45.0 


These values are plotted in Fig. 3-2. 

# Discussion 

The growth of the focused spot size with range is faster than the first 
power of the range and the growth of the focal depth with range is faster than 
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the range square, which will enable us to explain the signal-to-noise ratio 
at longer range falling off faster than l/(range) . 

3.3.3 Coherent Measurement 

• Previous Measurements 

The initial coherent spatial resolution measurements are plotted in Fig. 
3-3, When using SAT detector 952, the signal-to-noise ratio measurements 
at different ranges of that detector were not reliable at the moment the data 
were taken. As a result, the signal-to-noise ratio fell off unejfpectedly fast 
with range which is shown in Fig. 3-4. 

• Recent Measurements 

The signal-to-noise ratio as a function of range has been checked by- 
using Pb-Sn-Te detector 314 and another SAT detector PV 1031. The re- 
sults are shown in Fig. 3-5. The insufficient L.O. power and the possible 
lower quantum efficiency of the Pb-Sn-Te detector are responsible for the 
lower signal-to-noise ratio. But the behavior of the signal-to-noise ratio 
as a function of range should be preserved. In general, a ^2 dB variance on 
the signal-to-noise ratio is easily introduced into the results because of the 
laser power fluctuation, beam wandering, optical alignment, and the detector 
operation conditions. 

Signal-to-noise ratio is almost independent of the kind of sand target, 
i.e,, sandpaper wheel or belt sander. 

The coherent spatial resolution, as measured in the laboratory, is con- 
cluded in Fig. 3-6. 

3.3.4 Misalignment Effects on Spatial Resolution 

A Pb-Sn-Te detector was used in this test which had a 1.5-inch focal 
length Ge meniscus lens placed in front. The' range was 188 feet. 
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Fig. 3-3 - LiDV Coherent Spatial Resolution Measurements 
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Fig. 3-4 - LDV Signal- to- Noise Ratio Measurements as a Function of Range 


LMSC-HREC TR D390 159-1 



3“'5 ' LfDV Signal^to-Noise Ratio Measurements as a Function, of Range 


LMSC-HREC TR D390159-I 


Coherent Spatial Resolution (ft) 


I.MSC-IIREC TR D390159-I 



Fig. 3-6 - LDV Coherent Spatial Resolution Laboratory Measurements 
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• Defocused Detector 

The measured coherent spatial resolution was 10.8 feet operated under 
the backscatter mode. The range resolution was reduced to 9.4 feet when the 
detector was defocused by 0.20 inch either toward or away from the lens but 
was within the experimental error. 

# Lateral Displacement of Secondary Mirror 

The coherent spatial resolution increased from 11.5 to 12.5 feet when 
the secondary mirror was moved from the optimum position laterally 
toward the table center by 0.02 inch. The reflected transmission power from 
the primary mirror was shifted to one side of the optical axis when the above 
adjustment was made. 

# Misalignment of the Secondary Mirror Orientation 

Misalignment of the secondary mirror orientation did not affect the 
spatial resolution measurement by any significant amount when diffraction 
rings, which indicated good orientation of the secondary mirror, disappeared. 

• Misalignment of the Recombining Beamsplitter 

By moving the target, the spatial resolution under the separate L. O. 

mode was 10.5 feet at the optimum position of the recombining beamsplitter. 

The spatial resolution remained the same when the recombining beamsplitter 

-3 

was rotated by 0.81 x 10 radian and consequently lowered the signal 
level by a factor of 0.43, 
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3.3.5 Interference Effect Between the Backscatter Power and the Separate 
Local Oscillator Power 

# Experiments 

A potential problem with operating the LDV system in the separate 
L.O. mode is the interference between the backscatter power and the separate 
L.O. power which will give a fluctuating resultant L.O. power on the detector. 

The amount of fluctuation due to the scanning secondary — even when the ratio 
between the separate L.O. and backscattered L.O. powers is increased to 36 
to 1 — would still be 3 dB. 

The wheel signal level varies with the orientation of the recombining 
beamsplitter as indicated in Fig. 3 -7a. The translation of the recombining 
beamsplitter along the optical path (x'-axis in Fig. 3-8) was approximately a 
quarter wavelength. To vary signal level from maxima to minima, or vice 
versa, at each orientation of the recombining beamsplitter, the amount of the 
fluctuation, which is the difference between maximum and minimum, at each 
orientation is indicated in Fig. 3-7b. 

By translating the secondary mirror, the same type curves were obtained 
as shown in Figs. 3 -9a and 3- 9b. One-quarter wavelength travel of the secondary 
mirror changes the signal level from maximum to minimum. Rotation of the 

recombining beamsplitter about the z-axis (Fig. 3-8) by an amount of approxi- 

-3 

mately ^0.5 x 10 radian gives similar results which are plotted in Figs. 3- 10a 
and3-10b, respectively. 

® Discussion 

The ratio of the separate L.O. power level to the backscatter power 
level should be increased to at least 36 to 1 in order to cure the interference 
problem. 
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Wheel Signal Level at 188 Feet (dBm) 


SAT Detector PV 952 
Translation of Recombining 
Beamsplitter 


Note: Recombining beamsplitter position 

about Z'-axis = 0 rad. 


Recombining Beamsplitter Orientation 
(rotational adjustment about y-axis in mm) 


Fig. 3-7a - Interference Effect in the Separate Local Oscillator Mode 
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Fluctuation Amplitude (dB) 


Note: Recombining beamsplitter position 
about z-axis = 0 rad. 



Fig. 3 -7b - Interference Effect in the Separate Local Oscillator Mode 
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Fig.3“8 - LDV Interferometer Set-Up for Separate LiOcal Oscillator Operation Mode 
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Recombining Beamsplitter Orientation 
(rotational adjustment about y-axis in mm) 


Fig. 3“9a ” Interference Effect in the Separate Local Oscillator Mode 




Note: Recombining beamsplitter position 
about z-axis = 0 rad. 


, 1 .'ll 

SAT Detector PV 952 
Translation of Secondary 
Mirror 






Fluctuation Amplitude (dB) 


Note; Recombining beamsplitter position 
about z-axis = +0.5 x 10'^ rad. 



2.38 2.42 2.46 2.50 2.54 2.58 2.62 

Recombining Beamsplitter Orientation 
(rotational adjustment about y-axis in mm) 


Fig. 3-lOb - Interference Effect in the Separate LtOcal Oscillator Mode 
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The low frequency optical beat of the scanning secondary mirror sepa- 
rate L.O. mode, can be cured by the following arrangement: (1) introduce the 
separate ]L,0. onto a scanning flat mirror which is synchronized with the 
scanning secondary mirror; and (2) after the reflection by the scanning flat 
mirror, the separate L.O. can be introduced onto the recombining beamsplitter 
mixed with the return signal and the backscatter power? A frequency translated 
separate L.O. mode will not introduce the interference problem because there 
is no phase interference between the secondary mirror backscatter power and 
the frequency translated power when two frequencies are not equal. 

• Extended Measurements 

Note in Figs. 3-7a, 3-7b, 3-9a, 3-9b, 3-iOa, 3-10b, 3-1 la and 3-1 lb, 
that the maximum fluctuations of the wheel signal level did not occur when 
the wheel signal level was maximum and also that the top part of the wheel 
signal level was somewhat flat. One possible explanation of these results is 
the saturation of the spectrum analyzer because of the strong signal level. 

Switching the input attenuation of the spectrum analyzer from 10 dB to 
30 dB corrects the spectrum analyzer saturation. The result of the interference 
effect in the separate local oscillator mode is shown in Fig. 3-12 where it is 
seen that the top of the wheel signal is no longer as flat as before and the peak 
signal fluctuation occurs where the wheel signal is maximum. 

A 21 dB fluctuation in signal level was noticed. The destructive inter- 
ference gives lower resultant local oscillator power, hence the higher voltage 
across the detector and higher noise level, but higher signal- to- noise ratio 
that is shown in Fig. 3- 13. The constructive interference gives higher resultant 
local oscillator power, hence lower noise level but also lower signal- to- noise 
ratio, which is shown in Fig. 3-14. 




This technique is probably too complex to be practical. 
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Note; Recombining Beamsplitter Position 
about z-Axis = -0»5 x 10 rad* 


Recombining Beamsplitter Orientation 
(rotational adjustment about y-axis in mm) 


Fig. 3- 11a - Interference Effect in the Separate Local Oscillator Mode 


LMSC-HREC TR D390159-I 





Fluctuation Amplitude (dB) 


Note: Recombining beamsplitter position 
about z-axis =-0.5 x lO"^ rad. 



Fig. 3“ lib - Interference Effect in the Separate Local Oscillator Mode 
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Fluctuation Amplitude (dB) 



Fig. 3- 12b - Interference Effect in the Separate Local Oscillator Mode 

After the Spectrum Analyzer Saturation Has Been Corrected 
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Fig, 3- 13a - Sandpaper Wheel Signal Display of the Destructive Interference Between 
the Separate Local Oscillator Power euid the Backscatter Power (Range: 
188 ft; detector: SAT PV 952; Signal: +5 dBm; Noise: -71 dBm; S/N: 76 dB; 
Bandwidth: 30 kHz; Center Frequency 800 kHz; 0.1 MHz/division) 
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Fig. 3- 13b - Sandpaper Wheel Signal Display of the Destructive Interference Between 
the Separate Local Oscillator Power and the Backscatter Power (Range: 
188 ft; detector: SAT PV 952; Signal: +5 dBm; Noise: -71 dBm; S/N: 76 dB; 
Bandwidth: 30 kHz; Center Frequency 800 kHz; 0.1 MHz/division) 
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Fig. 3- 13c - Sandpaper Wheel Signal Display of the Destructive Interference Between 
the Separate Local Oscillator Power and the Backscatter Power (Range: 
188 ft; detector: SAT PV 952; Signal: +5 dBm; Noise: -71 dBm; S/N: 76 dB; 
Bandwidth: 30 kHz; Center Frequency 800 kHz; 0.1 MHz/division) 
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Fig. 3- 14a - Sandpaper Wheel Signal Display of the Constructive Interference Between 
the Separate Local Oscillator Power and the Backscatter Power (Range: 
188 ft; detector: SAT PV 952; Signal: -16 dBm; Noise: -76 dBm, S/N: 

60 dB; Bandwidth: 30 kHz; 0.1 MHz/division) 
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Fig. 3- 14b - Sandpaper Wheel Signal Display of the Constructive Interference Between 
the Separate Local Oscillator Power and the Backscatter Power {Range: 

188 ft; detector: SAT PV 952; Signal: -16 dBm; Noise: -76 dBm, S/N:60dB; 

Bandwidth: 30 kHz; 0.1 MHz /division) 
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Fig. 3- 14c - Sandpaper Wheel Signal Display of the Constructive Interference Between 
the Separate Local Oscillator Power and the Backscatter Power (Range; 

188 ft; detector: SAT PV 952; Signal: -16 dBm; Noise: -76 dBm, S/N: 60dB; 

Bandwidth; 30 kHz; 0.1 MHz/division) 
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• Signal- to- Noise Ratio vs Range for SAT Detector PV952 

After the saturation problem for the spectrum analyzer had been cor- 
rected, we found that either the local oscillator power was too high or the 
amplifier was not designed properly to provide optimized signal- to- noise 
ratio for that local oscillator power level which was accepted by the detector. 
This was noticed when the secondary mirror was moved later«d.ly approxi- 
mately 0.02 in. from its optimum position and the higher signal-to-noise 
ratio of the backscatter mode was measured because of the reduced back- 
scatter power from the secondary mirror due to this misalignment. 

At a range of 188 feet the highest signal-to-noise ratio from a sand- 
paper wheel was recorded as 80 dB. Under the normal operation, the signal- 
to-noise ratio was only about 75 dB for a total backscatter power of about 0.7 
mW accepted by the detector according to the I-V curve of each detector. The 
argument that lower local oscillator power gives higher signal-to-noise ratio 
was also verified by looking into the interference effect of the separate local 
oscillator mode. Destructive interference gives lower local oscillator power, 
which gives higher signal-to-noise ratio. 

The signal-to-noise ratio versus range for SAT detector PV 952 is shown 
in Fig. 3-15. In the figure the signal-to-noise ratio at a range of 188 feet was 
78 dB. This value was also reached by SAT detector PV 1031 once before, 
which encouraged us to further investigate the signal-to-noise ratio and shoot 
for higher values. 

3.4 FIELD MEASUREMENTS OF COHERENT SPATIAL RESOLUTION 

AND SIGNAL-TO-NOISE RATIO 

Spatial resolutions of the optical sub- systems of the laser Doppler 
units mounted in MSFC's vans 1 and 2 were measured at the field test site 
for comparison with theory and laboratory measurements during the month 
of June 1974. The results of these measurements are presented in Fig. 3-16. 
Signal-to-noise ratio (S/n) measurements were also made. The results of 
these tests are depicted in Figs. 3-17 and 3-18, 
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Fig. 3-16- Coherent Spatial Resolution Measurements versus Theory 
for LDVs 1 and 2 
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Coherent Signal-to-^Noi se Ratio, S/n (dB) Coherent Signal -to -Noise Ratio, S/N (dB) 


[>M.SC-HREC TR D 390 159-1 


VAN 1 



Fig. 3-17- Coherent Signal-to-Noise Ratio versus Target Range for LDV 1 


VAN 2 



Fig. 3-18 -Coherent Signal-to-Noise Ratio versus Target Range for LDV 2 
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Calculations of the signal- to-noise (S/n) ratio for untruncated Gaussian 
theory were made to compare with the actual system efficiency S/N ratio 
measurements. The results of these calculations and the comparison is de- 
picted as follows. 


Untruncated Gaussian Theory 

= Kap:i0:/R)2 


P • COS0- 
o 1 

hi/Af 
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Af 
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optics efficiency 
target reflectivity 
telescope radius 
range 

laser output power 

incident normal angle on target 

bandwidth 

photon energy 

Airy receiver efficiency 



Quanity 


System Number 


a = 


952, Lab 
0.24 

972, Van 1 
0.19 

1031, Van 2 
0.34 

p = 


0.07 

0.07 

0.07 

1 ^ ~ 


0.72 

0.30 

0.47 

r = 


0.5 ft 

0.5 ft 

0.5 ft 

R = 


200 ft 

200 ft 

200 ft 

P 
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20 W 

15 W 

20 W 

0 . = 
a 


45 deg 

45 deg 

45 deg 

Af = 


30 X 10^ Hz 

30x 10^ Hz 

30 X 10^ Hz 

hi' = 


1-87 X 10"^® 

1.87x 10‘^° 

1.87 X 10'^^ J 

K = 


-6 dB 

-6 dB 

-6 dB 

(s/n) 

power (calculated) = 

86.8 dB 

80.7 dB 

86.5 dB 

L . 

(measured) = 

80 dB 

75 dB 

75 dB 
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3.5 TESTS OF INTERFEROMETER AND CO^ LASER POLARIZATION 
3.5.1 Interferometer Alignment Tests 

Lockheed" Hunts ville personnel assisted NASA in evaluating the align- 
ment of the Raytheon fabricated Mach Zehnder interferometer {shown in the 
following sketch). The table below indicates the measured results of this 
test. These results imply that the system was well aligned when using the 
secondary backscatter as the local oscillator source, but was non-optimum 
when using the Mach-Zehnder configuration. 



Interferometer 
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Opei'ation Conditions 
200 Range, Bias Voltage = 0 

Bias 

Current 

Signal 

Noise 

s/n 


(mA) 

P-P(V) 

(dBm) 

(dBm) 

(dB) 

W/O Wire, W/O Separate L. O. 

5.0 

0.8 

+2 

-91.5 

93.5 

W/O Wire, W/O Separate L.O. 

3.8 

0.6 

0 

-91.5 

91.5 

W/O Wire, w/Separate L.O. 

7.0 

0.5 

-2 

-91.0 

89.0 

w/wire, W/O Separate L.O. 

0.5 

0.2 

-10 

-93.5 

83.5 

W/Wire, W/Separate L. O. 

2.5 

0.2 

-10 

-92.5 

82.5 


3.5.2 The Degree of Polarization Measurement of a Raytheon CO- Laser 
(Model SL-lOA) 2 

The following is a schematic of the experimental arrangement used to 
evaluate the degree of polarization of Raytheon CO2 laser. 


Coherent Radiation 
Model 201 Power Meter 
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Since the laser power was low (<^8 watts, deliberately adjusted), the tem- 
perature of the first Ge coupler was low. The absorption of the Ge coupler 
was therefore also low. Consequently, the laser power received by the 
coherent radiation power meter was approximately the same as that incident 
upon the first Ge coupler. The second Ge coupler ensured that the trans- 
mitted laser beam through that coupler was a p-wave (which appeared 
to the first Ge coupler as an s-wave). After the reflection by a high reflec- 
tance mirror, this s-wave laser beam was introduced into the model Cl 
thermopile (Sensors Inc.) through a calcium flouride attenuator (-12.7 dB) 
and a BaF^ lens. The output dc voltage of the thermopile was linearly pro- 
portional to the power of the s-wave laser beam for this lower power case. 
The calibrated responsivity of the thermopile enabled us to calculate the 
power of the s-wave laser beam from the first Ge coupler. The degree of 
polarization of the laser beam can be calculated as 


Degree of Polarization 


Total s-Wave Power 
Total Laser Power * 


where the s-wave power is much less than the total laser power. 

According to the experimental results the degree of polarization of 
the Raytheon CO^ laser was 99.20%. Another method was also employed 
to determine the degree of polarization. This method was a direct measure- 
ment of the s-wave power with a Quantronix Model 502 power receiver and 
a model 504 energy-power meter. This method yielded a value of 99.14% 
which was comparable to the value measured by the first method. 

3.6 DESIGN. CRITERIA FOR THE PHOTO VOLTAIC DETECTOR BIAS 

CIRCUITRY FOR CO^ COHERENT DETECTOR 

The design criteria of the bias circuit for the photovoltaic detector 
are outlined as shown on the following page. 
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• Safe operation v/ith the detector 

• Maximum photo-induced electrical power coupling 
from the detector to the preamplifier 

• Low noise, flat response preamplifier and circuit 

• Sufficient frequency bandwidth 

• Minimum radio frequency interference 

• Convenient operation 

Each design criterion is explained in detail as follows: 

3.6.1 Safe Operation With the Detector 


• Detector must be cooled using liquid nitrogen before 
the bias is applied. 

• Independent of power failure, optical local oscillator 
power blockages, or loose wire connections from the 
detector to the bias supply, the bias voltage generated 
across the detector must always be less than the max- 
imum allowable voltage applied to the detector. 

• The bias current must also be controlled by limiting 
the amount of optical power incident upon the detector 
when the safe bias voltage is applied. 

• A front panel light bulb is required on the bias box as 
a reminder to the operator to turn off the bias after 
each operation. 

• A means of achieving these safe operations is to use 
a small source voltage and small load resistance 
across the detector as illustrated below. 
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3.6.2 Maximum Photo-Induced Electrical Power Coupling from the Detector 
to I he Preamplifier 

• The equivalent circuit of the detector is 



Z 2 = R = series resistance (dc resistance of the 
detector) 

= load impedance 

Gj^ = small signal shunt conductance (slope of I-V curve) 
Cp = shunt capacitance 
CO = IF frequency 
i = detector generated IF current 

The above circuit can be redrawn as 
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where i' is called short-circuit current. In this circuit 
the load impedance will sense the same source current 
and impedance as in the previous circuit. 

• The condition of maximum power can be coupled by 

from the detector when Z 2 «Zj if Zj«z^. An amplifier 
having an input impedance equal to the detector impedance 
is preferred in order to provide maximum power coupling; 
however, because of amplifier availability, this is very 
unlikely. An alternative way to achieve the same purpose 
is to use a wide band transformer for matching the detector 
impedance to the input impedance of the amplifier. This 
circuitry is shown below; 


N;1 




The turn ratio, N, should equal to ^/z, for maximum 
power coupling. 


3.6.3 Low Noise Flat Response Preamplifier and Circuit 

The noise spectrum and the gain of preamplifier must be flat over the 
required frequency bandwidth for easier signal processing and calibration of 
the system. A low noise preamplifier and circuit is required for shot noise 
limited operation of the detector with a moderate optical L.O. power. 

3.6.4 Sufficient Frequency Bandwidth 

Both the preamplifier and the bias circuit must provide sufficient fre- 
quency bandwidth. Electrical component size, arrangement, and connections 
should be arranged carefully to avoid stray capacitance and coupling losses 
from the detector to the preamplifier. 
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3.6.5 Minimum Radio Frequency Interference 

RF shielding of the detector, bias circuit, and preamplifier must be 
included in the design of the detector bias circuit. This includes precaution 
against ground loops, RF isolation of the bias voltmeter and ammeter, and 
minimum wire lengths. For example, the bias voltmeter can be arranged 
as depicted below to avoid RF pick-up from the voltmeter: 


I 1 



3.6.6 Convenient Operation 

The bias box should be compact with bias voltage and current meters 
attached. The assembly should be designed for easy opticad alignment moni 
toring, optical power monitoring, and ready location of detector optimum 
operation conditions. Easy adjustment and fine control of the bias should 
be available. Convenient input and output connectors should be identified 
with labels. The box should be designed with no exterior exposed voltage 
terminals. 
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3.7 ACOUSTIC-OPTICAL TRANSLATOR EVALUATION 
3.7.1 Basic Principles 

When a laser beam, characterized by its wave vector, Ic, and its angular 
frequency, w., impinges on a traveling acoustic wave, w^, established in 
an electro-optic medium, the two waves interact to generate a third wave, kj . 
This interaction can be represented by a formulation of the conservation 
of momentum and energy 


kj = k. +k 
—1 a 


Wj = w, + w 
d 1 — a 


Figure 3-19 shows schematic and vector diagram of the interaction. The de- 
flected beam exhibits a frequency shift equal in magnitude to the frequency of 
the acoustic wave. 


This interaction, known as the Bragg effect, shifts the incident beam 
into two beams angularly separated by: 

A0 = Af/nv 

inside the elasto-optic medium, where 


A = the optical wavelength of the incident beam 
f = the acoustic frequency applied into the cell 
n = the index of refraction of the elasto-optic medium 
V = the acoustic velocity in the elasto-optic medium 
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(b) Down-Shifted in Frequency (c) Up-Shifted in Frequency 


Fig, 3-19 - Acousto-Optic Interaction 
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As shown in Fig. 3- 19a the deflected beam will be down- shifted in fre- 
quency. This down-shift is shown in Fig.3-l9b. Any amount of reflection of 
the sound wave, which lies in the same plaine as the incident sound wave, by 
the acoustic sink device will deflect the incident laser beam into the same 
direction as the beam which is deflected by the transmitted sound wave but 
up-shifted in frequency (Fig. 3- 19c). 

An optimum interaction is obtained when the angle of incidence of the 
laser beam onto the acoustic wave front is set at the Bragg angle 

q - M - AL 

B 2 2nv 

which is equal to 0.36 deg for X = 10.6 jtxm, f = 24 MHz, n = 4, and v = 5.09 x 10^ 
m/sec. 

Although the deflection angle, 9 deg, outside the Ge crystal is large because 
of the high index of refraction of the Ge, a precise adjustment of the Ge 
crystal orientation is still required for the optimum operation of the cell. 

3.7.2 Acoustical CO^ Laser Frequency Translator Test 

• Incoherent Measurement 

The experimental set-up is shown in Fig. 3-20. The frequency translated 
separate local oscillator power which is collected by the lens into the detector 
is shown in Fig. 3-21 for various knob settings of RF driver for the frequency 
translator . 

The efficiency of the frequency translator is defined as the ratio of de- 
flected laser power and incident laser power (see Fig. 3-22). 
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• Coherent Measurement 

The frequency spectrum of the sandpaper wheel signal Is shown in Fig. 

3-23. The 24 MHz carrier was a beat signal between the translated L. O. laser 
power and the untranslated laser power. Most of the untranslated laser power 
was the residual backscatter power from the telescope secondary mirror even 
when the attenuation wire was in place. As the secondary mirror backscatter 
was blocked, the 24 MHz beat signal dropped by approximately 42 dB. 

The 48 MHz beat signal was due to the beat between the laser beam which 
is up-shifted in frequency caused by the reflected acoustical wave from the 
acoustical sink and the laser beam which was down-shifted in frequency pro- 
duced by the transmitted acoustical wave. 

As the laser beam, focused by the telescope, was redirected from the 
lower part of the rotating sandpaper wheel to the upper part of the rotating 
wheel, a {24 - 1.4) MHz wheel signal and a (24 + 1.4) MHz wheel signal were 
seen as photographed in Fig. 3-24a and Fig. 3-24b, respectively. 

The "image signal," which is smaller in amplitude, is shown in Figs, 3-24a 
or 3-24b was a result of the wheel signal beat with the L.O. power up-shifted 
in frequency which was produced by the reflected acoustical wave from the 
acoustical sink. The signal strength difference between the "signal" and the 
"image signal" can be adjusted to as high as 30 dB with extra alignment care. 

The noise spectra of the backscatter mode and the frequency treuislated 
mode are shown in Figs.3-25a and 3-25b, respectively, with the telescope output 
mirror blocked. 

The lower noise level and the spikes shown in Fig. 3- 25b is due to the strong 
24 MHz signal saturating the detector preamplifier and over-driving the spectrum 
analyzer when the input attenuation of the spectrum analyzer was 0 dB. The 
noise spectrum is cleaner with higher input attenuation of the spectrum analyzer 
which is shown in Fig. 3-26. 
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(a) Telescope Output Mirror 
Unblocked 



0 Hz 




100 MHz 



(b) 


0 Hz 


100 MHz 


Telescope Output Mirror 
Blocked 


Bandwidth: 30 kHz 

Spectrum Analyzer Input Attenuation: 20 dB 
LOG REF Level: 20 dBm 

Range of Target: 71 ft (Folding Mirror is Used) 

Fig. 3-23 - Frequency Spectrum (0 to 100 MHz) of the Sandpaper Wheel Signal 
with Frequency Translated Separate Local Oscillator Power 
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24 MHz 

(a) Scattering Particles of the Sandpaper Wheel Move Away from the 
Incident Laser Beam 



-H 1 MHz h- 


24 MHz 

(b) Scattering Particles of the Sandpaper Wheel Move Toward the Incident 
Laser Beam 


Fig. 3-24- Sandpaper Wheel Signal in the Frequency Translated 
Separate L.O. Mode 
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«*- _40 flUm 


0 Hz 


50 MHz 


Input Attenuation: 10 dB 


(a) Backscatter Mode 



(b) Frequency Translated Mode 


Fig. 3-25 - Noise Spectrum with Telescope Output Mirror Blocked 


Receding page blank not cilmeo 
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Fig. 3-26- Noise Spectrum of the Frequency -Translated Separate L.O. 

Mode with Telescope Output Mirror Blocked and Higher Input 
Attenuation of the Spectrum Analyzer 
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The backscatter mode (untranslated frequency mode) wheel signal and 
noise are shown in Fig. 3-27. 

The frequency translated mode wheel signal and noise are shown in 
Fig. 3-28. The signal-to-noise ratio measurement in this mode was not valid 
because of the preamplifier saturation caused by the extremely strong 24 MHz 
beat signal between the frequency translated L. O. beam and the residual 
backs cattor power from the secondary mirror. 

3,7.3 Suggestions 


• Cut the back surface of the Ge acoustical modulator crystal, 
which is opposite and parallel to the surface with transducer 
attached, into a wedged surface to eliminate the collinear 
reflection of the acoustical wave from that surface. 

• Install a low-noise notch filter with enough attenuation and 
matched impedance in front of the preamplifier to overcome 
the amplifier saturation by the 24 MHz strong beat signal. 

• Install an adjustable mount under the acoustical frequency 
translator. 

• Install a lower RF power driver for the photovoltaic detector. 
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(a) Signal 


1 MHz 

f 

1.2 MHz 


+ 10 dBm 


Input Attenuation: 30 dB 
Bandwidth: 30 kHz 



(b) Noise 


f 

1.2 MHz 


30 dBm 


Input Attenuation: 10 dB 
Bandwidth: 30 kHz 


Fig. 3-27 -Backscatter Mode Sandpaper Wheel Signal and Noise Display 


161 



r.M.SC-]fH]*X; TK D390159-I 



1 MHz 

(24-1.2) MHz 


(a) Signal 


+ 10 dBm 


Input Attenuation: 30 dB 
Bandwidth: 30 kHz 



(24-1.2) MHz 


(b) Noise 


— -50 dBm 


Input Attenuation: 10 dB 
Bandwidth; 30 kHz 


Fig. 3-28- Sandpaper Wheel Signal in Frequency Translated Separate L. O. 
Mode as well as the Noise Display 
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Section 4 

LASER DOPPLER VORTEX SENSOR SIMULATION PROGRAM 

Authors: K. Shrider 
E. Coffey 


4.1 BRIEF DESCRIPTION 

The program developed simulates the output of NASA's laser Doppler 
aircraft wake vortex sensor. 

Simulated are the following: 

0 Air velocity due to the wind and vortices 
o Vortex transport 

© Particle size distribution in the atmosphere 
© Backscatter from particles in the illuminated volume 
o Sensing volume motion due to scanning 
® Detector response 
© Filter response 
o Filter integration time 

& Digital tape output and printout similar to that of the 
real data processor. 

4.2 DETAILED DESCRIPTION 

The simulated parameters are now discussed in more detail. 

• Air Velocity 

The velocity field in the scanned plane is composed of two vortices and 
their ground images, superimposed on the wind field. The plane is perpendicular 
to the axes of the vortices. No wind turbulence model is currently anticipated 
for this program. 


PAGE BLANK NOT mMm 
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The vortex circulation, is computed from the aircraft's speed, 

weight, and wing span as follows: 


ro 


4w 

ttP u s 


(4.1) 


where p = density of air, w = aircraft weight, u = aircraft speed and s = wing 
span. 


The velocity field of a vortex is represented by 



where 

r = the distance from the center of the vortex 
t = the time from vortex creation 
U = the eddy viscosity. 

The wind field is a curve fit of a set of velocity and direction data points 
as a function of altitude. 

• Particle Distribution 

Two types of particle species can be accommodated and are treated 
separately. Both size distributions are of the following basic shape (see 
also sketch below) at a sufficient distance from the vortex core. Near the 
core the centrifuge effect is accounted for. 
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where 


(a) 

N(a) 


N 

C 


N(a) 



a j < a < a^ 

a 2 < a < a^ 
a 3 < a 


the radius of the particle in micrometers 

the number of particles per cubic centimeter 
within a 1 mm radius increment centered about 
a radius a 

the total number of particles per cubic centimeter 
a normalizing factor equal to 



(4.3) 


The redistribution due to the centrifuge effect is a function of age, the velocity 
distribution of the vortex, the radius and density of the particles, and the tem- 
perature and density of the air, 

• Backscatter Cross Section, P 

P , the equivalent backscatter cross section of the particles per unit 
volume, is computed from the Mie series and is a function of laser wave 
length, particle size, particle index of refraction, and the particle number 
density. 

• Detector Response 

The rcinge dependence of the response of the LDV is given by 


tan (L 2 - L^/Lj^) - tan (^1 " 
ir/2 +tan’^ (L^/Lj^) 


167 



LMSC-HREC TR D390159-I 


Line of Sight 


Transceiver 


Range 


The integration process along the line of sight of the LDV is carried out by- 
summing over length elements Lj L 2 where 


L 


the fraction of the detected power at the LDV output 
arising from the volume between a plane at Lj^ and 
a plane at L^ 

range of maximum response (approximately the focal 
range, f) 




1 + 


(XfVjra^) 


where X = the wave length of the laser radiation 

a. = the exp(-2) intensity radius of the laser 
beam at the primary mirror. 


The above response function holds true for an infinite mirror without 
any obstruction. To correct for the finite mirror size and the obstruction 
in front of the mirror the quantity a^ is substituted for aj^, thus, 

z z 

^e ” ^b ^ 


where 

= radius of obstruction due to the secondary mirror 
= radius of primary mirror* 


The contribution to the detector output power from between Lj and is: 


^(Lj.L2) ^ ^Lj,L2) ^ 


(4.6) 
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where 


{Li, L^) 


P 


the fraction of received power from the volume 
between and 

the fraction of the energy backscattered by the 
particles per unit volume 

a normalizing factor given by 




(4.7) 


to ensure that P 




has units of power . 


where 


n 


X 

h 


c 


the quantum efficiency of the detector 
the output power of the laser 
the wavelength of radiation 
Planck's constant 
the speed of light 


# Filter Bank Response 

The illuminated volume is partitioned into contiguous segments (volume 

between and such that the velocity of all the particles In a particular 

segment can be considered to be traveling at the same velocity. The power 

P T i \ from the segment is represented by a Rayleigh distributed random 
(LiFf-'i+f) 

variable to account for phase cancellation caused by the arbitrary separation 
of the particles in the segment. 


The line of sight velocity is computed and ccmverted to a Doppler shift. 

This frequency is compared to the center frequencies of the filters to locate 

the signal on the characteristic curve for each filter. The magnitude of the 

response from each filter is computed by multiplying P^^^ j by the value 

1 2 

of the characteristic curve at the Doppler shift frequency. Also, for each 
filter response, a random phase angle is assigned. 
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For each filter the response vectors caused by the illuminated volume 
segments are added together. Also added is the random shot noise. 

The magnitude of the shot noise in each filter is a Rayleigh distribution 
of the RMS noise where 


RMS . =Jn \B/hc 
noise \ ^ 

where 

^ = detector qaantum efficiency 

X = filter bandwidth 
h = Plank constant 
c = velocity of light 
B = filter bandwidth. 

Also a random phase shift is assigned. 


(4.8) 


• Signal Integrator 

Each filter is sampled at 15.625 x 10 ^ second intervals. The integrator 
adds for each filter the sampled magnitude over a specified interval. This 
summation for each filter is multiplied by a gain factor where the gain factor 
is given by 


Gain 


8 

(RMS Noise) * (Number of Samples in Each Summation) 


(4.9) 


During the integration some correlation of filter output magnitudes is 
expected since some of the sampled returns are from the same set of particles. 
Therefore at each partition the cross beam particle velocity relative to the 
beam is computed. Then, from that computation, the transit times for the 
particles to cross the beam width are computed. During the transit time the 
magnitude of signal from a partition is held constant and the phase angle is 
rotated as a function of the line-of-sight velocity. 
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• Laser Beam Scanning 

The range scanning and angular scanning of the focal point are inde- 
pendent of each other. The focal point moves at a constant rate in range 
and moves at a constant rate in angle. 

Vortex Movement 

A vortex moves with a velocity determined by the sum of the velocities 
at its center due to its ground image, its companion vortex with its ground 
image, and the wind, 

4.3 FLOW OF PROGRAM 

Reference is made to Fig. 4-1. Major dots • indicate the instantaneous 
location within the program (i.e., the particular routine) and minor dots • 
indicate operations within that routine, 

• MAIN Program 

Call INPUT. 

• INPUT Routine 

• Reads input (see Section 4.4). 

• Computes equivalent beam radius at the primary mirror (Eq. (4,5)). 

• Computes preliminary wind calculations for subroutine WINDY. 

If desired, it calculates a "least square" curve fit of wind magnitude 
to a power law equation. 

WSPEED(I) = WSPR C POWER (4.10) 

where WSPEED(I) is the wind magnitude at altitude ALT(I). 

These are input along with ALTR (a reference height). 

WSPR and CPOWER are computed. 
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Fig. 4-1 - Flow Diagram of LDV — Vortex Interaction Simulation Program 
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Also, routine LSQPE is called to compute the least squares 
curve fit of the wind direction as a function of altitude. This 
curve is a polynomial. A polynomial least square curve fit 
is made to represent the temperature field. A pressure and 
density profile are computed from the pressure at the ground 
and the temperature profile. This is done by integration using 
routine RUNGK. 

• Vortex circulation is computed using Eq. (4.1), where the 
atmospheric density is obtained by interrogating routine ATM 
with the vortex height. 

RHO = atmospheric density obtained by interrogating routine 
ATM with the vortex height. 

• Eddy viscosity, UNU, is computed from the temperature, density 
and recirculation using an empirical relation obtained from the 
vortex decay literature, thus 


UNU = 


UUU 

RHO 


* 10. ** 


(■ 


62 * LOGjq ^GAMMA * 


^ RHQ i 

lUUUj 


•) 


1,06 (4.11) 


where UUU 


/ .3048 \ 1.458E-6 * (TH * S./9.) ** 1.5 

V14. 5939029/ [TH * 5/9 +ll0“74) 


TH is temperature in degrees Rankine from routine ATM. 

• Vortex separation is computed assuming elliptically loaded 
wings, thus 

Vortex separation = (jr/4) * WSPAN. (4.12) 

• Pj., a normalizing factor is computed using Eq. (4.7). 

• RMS noise is computed using Eq. (4.8). 

• A gain factor is computed for scaling the output using Eq. (4,9). 

• Other initializations and unit conversions are made. 

• Call function BATA. 


• BATA Function 

• Call SPHERE 
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SPHERE Routine 


• Computes the backscatter cross section for a single particle 
using the Mie series (a function of radiation wave-length 
particle size and index of refraction). 

• Return to BATA, 


BATA Function 


• Compute the backscatter cross section per unit volume for the 
particle distribution. Also, it stores backscatter cross sections 
for discrete size particles within the distribution. These are 
used later in conjunction with the centrifuge routine. 

• Return to INPUT. 


INPUT Routine 


• Routine RANDUM is called to initialize the random number 
s equence. 

• Call VORPOS. 

VORPOS Routine 


• Integrates for the vortex position. 

• Calls routine WINDY for wind velocity, 

• Calls for the integration routine RUNGK. 

• Uses standard technique for superimposing air velocity from 
the vortices, their ground images, and the wind. The velocity 
at the center of a vortex (neglecting that vortex) is assumed 

to be the transport velocity of the vortex. The tangential 
velocity due to a vortex at range, r, from the core is given 
byEq.(4.2) ® 

• Return to INPUT. 


INPUT Routine 


• Call SCAN. 
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• SCAN Routine 


• Computes the proper location for the focal point. The angle 
and range movement are independent. Both move at a constant 
rate. 

• Return to INPUT. 


• INPtrT Routine 

• Call OUTONE. 


• OUTONE Routine 


• Writes the first record on tape. 

The output unit is 12. 

The record is written in unformatted FORTRAN. 
The record contains 5 words. 

1. NRUN is the run number — Integer 

2. PRTIME is the integration time — Read ^ 

3. IVMODE is the velocity coverage mode — Integer 

4. LINEAR/CFAR Code - Integer* 

5. TLR CODE - Integer 

• Prints run number and integration time 

• Return to INPUT. 


• INPUT Routine 


• Return to MAIN. 


• MAIN Program 


• Determines whether any part of the focal volume is within a 
distance, VSEP (where VSEP is the initial separation) of the 
vortices. This is to isolate regions of potentially large velocity 
gradient (see PART below). 


Variables introduced by NASA-MSFC into their readout routines. 
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• Routine SCAN is called for focus position. 

• The location of the focal point is saved periodically (normally 
every .002 sec). At output time the last saved point is printed. 

• Call PART. 


• PART Routine 


• Calculates Gaussian beam parameters. is the range to the 

maximum intensity point; Lp is the distance along the axis 
from the point such that 1/2 the returned power is from the 
interval (L^ - + Lj^)- 

• Calculates a cell size along the axis of the beam centered at 
the Lq point. This cell will contain a predetermined (typically 
.2) fraction of the power collected by the detector. Then if the 
focal point is close to the vortices, VORVEL will be called to 
yield the air velocity at the end points of the cell. If the velocity 
gradient is too large the cell will be cut in half. Cuts will be 
made until the velocity gradient is not too large. 

• Calculates a centroid for the cell weighted by the power returned 

to the detector. In the case of the first cell the centroid is L . 

o 

• Calculates the fraction of returned power which comes from 
this cell. 

• Calculates the beam width at the centroid. 

• Call RETPOW. 


• RETPOW Routine 


• Call VORVEL. 


# VORVEL Routine 


• Call WINDY. 


• WINDY Routine 


• Calculates wind velocity at specified height 

• Return to VORVEL. 
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• VORVEL Routine 


• Computes air velocity at a specified point (centroid of cell) 
due to the two vortices and their ground images, then adds 
the wind to get the total air velocity at the point. 

• Return to RETPOW. 


• RETPOW 

• Computes the cross beam air velocity at a cell centroid 
relative to the beam (includes scanning motion) 

• From the cross beam air velocity and the beam diameter 
at the cell centroid the transit time for particles to cross 
the beam is computed. 

• Computes the line-of-sight air velocity at the cell centroid 

e From the line-of-sight velocity the Doppler shift is computed. 

• The range from the cell centroid to the nearest vortex is 
computed. (This is to account for the centrifuge e^ect.) 

• If the centrifuge effect is important call BATA. 

• Call function BATA. 


• BATA Function 


• If routine CENTRI is called within a certain interval since it 
was last called (specified in input) it is skipped. (The centri- 
fuge effects needs to be updated only at time intervals large 
in comparison with the cycle time of other computations.) 

• Call CENTRI or CENTR 

The entry point CENTRI is for the full centrifuge representa- 
tion and the entry point CENTR is for a short approximation to 
the centrifuge effect. 


• CENTRI Routine 


9 Computes the radial particle size distribution from the center 
of a vortex as a function of vortex circulation, temperature, 
density of air, density of particles, size of particles, and time 
from creation of vortex. 

9 In the case where entry point CENTR is used only the radius 
of the empty core for one particle size (predetermined in the 
input) is computed. 
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• Return to BATA. 


# BATA F unction 


• Using routine CENTRI the backscatter area per unit volume 
is computed at 250 radii from the center of the vortex. 

• The backscatter is then computed for the cell centroid using 
a polynomial curve fit of three of the 250 points discussed 
above. 

• In the case of the simplified centrifuge model the core is 
assumed to have no backscatter and outside the core is 
assumed to have the same backscatter as the undisturbed 
air. 

• Return to RETPOW. 


• RETPOW Routine 


• The power output from the detector originating from the cell 
is computed using Eq. (4.6). 

• Call FILTER. 


• FILTER Routine 


• A random detector voltage is selected from a Rayleigh distri- 
bution about the mean power computed in routine RETPOW. 

The random voltage is to take account of power fluctuations due 
to phase cancellation due to the random motion of the particles. 

« Using the Doppler shift frequency and the characteristic curve 
of the filters the amount of voltage in each filter due to a cell 
is determined. 

• A random phase angle is applied to each filter output originating 
from the cell. The voltage is resolved into components for sum- 
mation over all the partition cells. 

• When this routine calls for the last cell, a random Rayleigh 
root -mean -square noise is added to each filter summation. 

This is to simulate shot noise. 

• Return to RETPOW. 
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• RETPOW Routine 


• Return to PART. 


• PART Routine 


• The next cell in the partitioning which extends from the last 
cell boundary is now accounted for and routine RETPOW is 
called. This process continues, the range being increased 
until the outer boundary of a cell exceeds a range equal to 
Eq + 3 Lj^. If the line-of-sight airspeed is increasing, cells 
will be constructed out to a rajige L.^ + 6 Lr. In no case will 
a cell extend beyond a range Lq + 6 Lj^. At this point cells 
are constructed inward from the first cell {toward the laser) 
to a range - 3 Lj^ or Lq - 6 

• Return to MAIN 


• MAIN Program 

• Time is increased. 

9 If sufficient time has elapsed, routine VORPOS is called to 
reposition the vortices. 

• Routine SCAN is called to re-position the focal point. 

• If the focal point has not moved far enough to justify re-, 
partitioning then routine FILTER is called directly, 

9 Call FILTER. 


9 FILTER Routine 


9 If the particles at the centroid of the cell have not departed 
the beam since the routine was last called, the phase angle 
is rotated by a function of the line-of-sight velocity. 

9 If the particles at the centroid of the cell have departed the 
beam, a new random voltage is selected from the Rayleigh 
distribution (about a mean power originating from the cell). 
Also, a new random phase angle is selected. 

9 The voltage in each filter due to the cell are summed and noise 
is added. 

9 Now the absolute value of the voltage in each filter for the 
current cycle is added to the time of the parameter for 
previous cycles in the integration time. 
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• Return to MAIN. 


• MAIN Program 

• If sufficient time has passed, routine OUTPUT is called. 

• Call OUTPUT. 

• OUTPUT 


• The gain is multiplied by the summed power in the filters. 

• Writes all but first record on tape unit 12. The records 
are written in unformatted FORTRAN. The logical records 
contain 107 words of 36 bits. 

Word 1 = frame number — Integer 

2 = x-coordinate of focal point — Integer 

3 = y-coordinate of focal point — Integer 
4-107 = contents of the 103 filters — Integers 

• A printout of this record is made. 

• The summed filter power is zeroed. 

• Return to MAIN. 


• MAIN Program 

• The time is increased and the program continues in the 
above manner. 
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4.4 VORTEX SCAN PROGRAM INPUT GUIDE 

TiTLt CARU IZ COLUiANS OF COIhMcNT 

(IF F IRbT 6 COLOMNO ARE dLANK» PROGRAi'-i TERMINATti) 


AiRGRAFT CARO CHOObt STANUARU AIRCRAFT ih I TH FIwST 6 COLOi''tNb 
(LEFT JUSTIFIED) 

ThL CHOlCtS are b/47 B7Q7 U7E7 DCIQ DCV DCS CVSoO 
LiOll CbA 

ace Et_ANtcS INDICATES AIRCRAFT FROr-i PREVIOUS CASE 
W I LL EE USED 


NR UN 
EPS 1 


LPSE 

PAR MAX 
NR AN 


naFiEL I ST 


( RUN ) 



Kui\| number 

CK I TER I UN FOR CEcL DIVISION 

FK ACT ION OF RETURN RAuIATION 

CK I TER I UN FOk cell D I V I S I ON 

DELTA VELOCITY IN CELL (METERSXSEC) 

MAXIMUM CELL SIZE (METERS) 





staktino Number for random number generator 


namelist (VORTEX) 


I VORTX 
Ru- IGHT 
SPEtD 

wE I GHT 
vVSPAN 
A I RHU 

xvort 


number of vortices (CURRENTLY ONLY CHOICE IS 2) 

VORTEX ALTITUDE ABOVE GROUND (FT) 

speed of aircraft kELATIVE to GkOuND (FT/SEC) 

IF NEGATIVE« (SPtED) lo TAKEi'I AS AIR SPEED 
WEIGHT OF AIRCRAFT (POUNDS) 

WING SPAN OF AIRCRAFT (FT) 

AIRCRAFT HEAuiNG IN DEGREES CLOCKW 1 St FROM 
magnetic NURTH 

CROSS Runway position for centek of vortex pair (m) 


Namelist (winds) 


JW I ND 


UW I NDT 
W 1 NOD 
i'lW 

<«LT { 30 ) 
WSPtED ( 30 ) 
WD IREC< 30 ) 


FLAG 

JWlND = 1 WIND IS PROVIDLDBY TAPE 10 AS A 

FUNCTION OF time 
(CURRENTLY NOT AVALAoLE) 

jWiND = 2 WIND IS PROViDEu oY TABLE wiTM tlNciAR 
interpolation AiNU linear EXTRAPi^tAT IOIn 
AS A function Uh altitude 

JWIND = 3 wind is PROVIuED BY LEAST SOUARe CURVES 

OF table values 

TIMt J.N WIND TAPE WHICH COkkESPONDS TO ZERO TIME OF 
SIMULA! I ON 

direct I UN OF wind reference ON TAPE IN DEGREES 
CLOCK-WISE FROM ( -Y ) 

NUMBER OF ALTITUDES 

altitudes at which the WINuS are specified (FT) 

WIND SPEED (FT/StC) 

dirlCtion f Rum which w I nu is Coming* 
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\UL,\J ( JO ) 
Ml- r>^ 

WSpK 

CPOv-jEk 


I'JPoly 
C on F < Sy ) 


PPEiiSG 

JTEMP 


|>I A 


aL.T it (JO) 
T!-i«1P t JO ) 

• JHDL YT 

(y ) 


Ot-GktLJ Akt Mc.AJUKt.LJ CLUUkwlJt FkUM MA uNc-T 1 C 
NUkTH 

VLKTICAL WINU JkcihU (KT/SLU) 

KLFt-kENCL AuTITUUt F Ok LEAST JoUAkU CUKVL UF 
SPEED (FT) 

iVliMD SPEED AT kEFEkENCE ALTITUDE (FT/SEL) 

< COMPUTED- WFIEk Ni*/ IS Gkt.ATEK THAN 0) 

Is VAklABLE IN LEAST SUUAkt^S CUkVE 
WSPEED/WSPP = ( ALT/ALTk)*^tCPUWEk 

(COMPUTED WHEN NW IS GkEATEk THmN 1 ) 

UKUEk uF POLYNOMIAL TO BE FlTTn-D TO WIND DlkcCTION* 
MAXlFiUF't VALUE FOk NPOLY IS 7 
Coefficients of polynomial fuk aIind jIREctiun 
(COMPUTED WHEN NW IS GKEATEK THAN NPOLY ) 

PKESSUkE AT OkOUND uEVEL IN MlLuIoAkS 
FLAG 

JTDMP = 1 TEMPEkATUkE IS PkOVlUED UY TABLE WITH 

LlNtAk IiMTEkPOLAT ION A(SU cINEAk 
EXT kAPOLAT ION AS A FUNCTION OF ALTITUDE 
JIEMP = 2 TEMPEkATUkE IS PkOVlDED oY A LLAbT 
SOUARES polynomial uF TAbLE VALUES 
NUi-ibLk '.jF A(_T 1 TUuES AT wH i CH THu TLMPEkATUKc. IS 

specified. 

ALfITUUES AT 1 ILH THE TEi'iPEKATUkE Is SPEC I F I eD . ( K T ) 
TLMPEkATUkE I DEGkEES C 

OKDLk OF PuLYinOi-'i I AE TO BE F ITTeu TO TEMPeRATURE 
MAXIivUm value Fuk NPOLY F IS 7 
Coefficients of polynofiIal fok lEi’-'PEKATUkE ppofiee 
(computed when na ib (jPeatek Than npolyT) 


namelist (LASEk) 


PowEk 

eJANDW 

wave 
kAjlBE 
KAJ I OT 
k 1 I N 
LDANEF 


LASEk OUTPUT Cw POwEft (WATTS) 
filter bANuwIdTH (HL) 

WA\/E length of LASEk kADlATION (McTEkS) 

KaDIUs of laser oEAM AT PkI|viARY LENS (mETERb) 
OJTEk RADIUS OF PKImAwY LENS (meTERS) 

I.MNER RADIUS OF PRIMARY eENS (MeTERS) 

DUANTUM efficiency ( D I MENS 1 ONLeSS ) 


namelist (SCANER) 


YL-ASeR 

X.-ASER 

FvJCALL 
^ ECALU 
THlTAL 


height of SOANEk FkOM GkOuNU (MeTEPS) 

CKoSS kUNWAY POST I ON OF SCANtk (METEkS) 

POSITIVE IS TO POkT OF AlkCkAFT 
MINIMUM kANGE F’Uk SCANEk (I'lETEkS) 
maximum kANGE FUk SCANEk (MlTEPS) 

MINIMUM ANGLE ABOVE HORIZONTAL FOR SCANLR (DEGREES) 


182 



LMSC-HREC TR D390159-I 


I HD Au 
t-ui u(- 

t- '■<!,. u r 
r»C AL 

sthlta 

Si GNri 
SiGHT 


NTyPE 
I ■'JDt- X ( 2 ) 

PWHO i Z ) 
Alt 2) 

fkZ ( 2 ) 

A j ( 2 ) 
Xi'JUMBK{2) 


Lit T A 


Pt,lZtiT 


ST I WL. 
i- T I AI2 
VST I ivit. 
SST ii'lE 
VsT I AiE 
pur IML 
SAiT K-iL 
Ul,T I ML 
,VPP 1 NT 


MAXIMUM ANGLE AbOVt- HUWl2UNTAb h UN SCANLN t LLGKlLS ) 
Fi^LUULNCY Ot NAIsiGt SLAN ( CYGLLt./SLL > 

FKtUUENCY UK angle scan (LYLLES/SLC) 

SI'AMTInG RANGE UF KUCAL POINT (METERS) 

starting angle from horizontal for focal point (DO) 

starting direction of range scan* 

+i«0 FOR POSITIVE -1.0 FOR NEGATIVE 

starting direction of angle SCmN. 


+i*0 FOR POSITIVE -1*0 FOR NEGAT 


namelist ( PARTI C) 


NUMBER OF PARTICLE TYPES (1 OR 2 ) 

I liiDEX OF refraction FOR EACH TYPE OF PARTICLE 
INDEX IS COMPLEX 
DENS ITT OF PARTICLES nG/K** 3 
LOvvER LIMIT OF PARTICLE RADII (M) 

CONSTANT D I -ilK I but 1 ui'* BET WEEN Al AND A2 
KeCIPRuCAE of FOURTH POWER BETWEEN A2 ANO 

Total number of particles 

FLAG 



AO 


JCENT 

= 0 

NO 

centrifuge 

NO 

ivi i E 

SCATER 

OCENT 

1 

NO 

centrifuge 

YLS 

Mid 

SCATER 

JCENT 

= e 

YES 

CENTR 1 FUGE 

NU 

Ml t. 

SCATER 

JCENT 

= 3 

YES 

CENTR IFuGE 

YLS 

Nl E 

SCATER 

JCENT 

= 4 

FULL 

CENTR IFUGE 

YES 

MIE 

SCaTER 


Eo J I VAEt.NT BACK SCATER AREA OF PARTICLES M**2/M**3 
BETA IS COMPUTED IF OCENT = 1 * 3 OR 4 

For centri f'oge. raoius in meters 
pakticle size used in determining step function 
(only used when OCENT = 2 OR O > 


NAMlLIST (TIME) 

PKOCciSSOR start time (SEC) 

PKOCESSOR CUTOFF TIME (SEC) 

VOPTEX START TIME (VSTIME MOOT BE *eT. ST I mE ) (PtC) 
SCAN START TIME (SSTIME MUST BE .EE. STIMD (SEC) 
STEP SIZE FOR VORTEX POSITIONING ( ^EC ) 

InTlRVAL between OUTPUTS (SEC) 
filter SAMPeING TIMl. 

InTbRVAL between position sampling for OUTPUT (SEC) 
NU.-iBER OF OUTPUT FRAMES TO EACH PRINTOUT 
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PRESET INPUT VALUES 


NAMELIST NAME VARIABLE VALUE 

run NRUN 

EPSl 0.2 

EPS2 0.5 

PARMAX 2.5 

NRAN 1.0 


VORTEX 

IVORTX 

HEIGHT 

SPEED 

WEIGHT 

WSPAN 

AIRHD 

XVORT 

WINDS 


JWIND 

2.0 

SWINDT 

0.0 

WINDD 

0.0 

NW 

2.0 

ALT{30) 

0., 10000.0 

WSPEED(30) 

0.0 

WDIREC(30) 

0.0 

WSV(30) 

0.0 

ALTR 

23.0 

WSPR 

0.0 

C POWER 

0.27 

NPOLY 

2.0 

COEF<9) 

0.0 

PRESSG 

1013.0 

JTEMP 

1.0 

NA 

2.0 

ALTIT(30) 

0.0, 2000.0 

TEMP(30) 

20.0 

NPOLYT 

1.0 

COEFT(9) 

0.0 


LASER 

POWER 

BANDW 

WAVE 

RADI BE 

RADIOT 

RADIIN 

QUANEF 


12.0 

120000,0 

0.0000106 

0.15 

0.15 

0.01905 

0.25 


2,0 
250.0 
-247.0 ) 

565000. 0> for B747 
195.7 ) 

0.0 
100. 0 
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SCANER 


PAR TIC 


TIME 


YLASER 

XLASER 

FOCALL 

FOCALU 

THETAL 

THETAU 

FREQF 

FREQT 

SFOCAL 

STHETA 

SIGNR 

SIGNT 


NTYPE 

INDEX (E) ( 

PRHO{2) 

Al(2) 

A2(2) 

A3(2) 

XNUMBR(2) 

JCENT 

BETA 

PSIZET 


STIME 

FTIME 

VSTIME 

SSTIME 

VDTIME 

PR TIME 

SMTIME 

DLTIME 

N PRINT 


0.0 

o.b 

40.0 

250.0 

5.0 

35.0 

6.0 
0.1 

40.0 

5.0 

1.0 
1.0 


2.0 

.212, 0.06), (1.650,0.09) 
1000,1192 
.02E-6, .02E-6 
.lE-6, .lE-6 
lOE-6, lOE-6 
1.378E+ 9, 1.378E+9 
0.0 

5.601492E-8 

1.35E-6 


0.0 

20.0 

-5.0 

0.0 

1.0 

0.002 

15.625E-6 

0.0002 

1 
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4.5 OUTPUT DESCRIPTION 

Two types of output are possible, a line printer record and digital mag- 
netic tape. 

A typical line printer output is shown in Fig. 4-2 and includes: 

• Integration time (milliseconds) 

• Time (seconds) 

• X (horizontal position in meters of focal point from LDV) 

• y (vertical position in meters of focal point from LDV) 

• RADIUS (range in meters to focal point from LDV) 

• THETAD (beam elevation angle — degrees ) 

• The column under CELL is the filter number. 

• The column under AMP is the integrated power from 
the corresponding filter. 

The digital magnetic tape is similar to that produced by the real data 
processor program. The tape contains record one, which has five words 
written in unformatted FORTRAN form. The five words are; 

1. Run number — Integer 

2. Integration time — Real 

3. Velocity coverage mode — Integer 

4. Linear/CFAR code — Integer 

* 

5. TLR code — Integer 

Record two, three, etc., contain 107 words written in unformatted FORTRAN 
form. The 107 words are; 

1. Frame number — Integer 

2. x-coordinate of focal point — Integer 

3. y-coordinate of focal point — Integer 

4. 107. Integrated voltage from the 103 filters. 

Variables introduced by NASA-MSFC into their readout routines. 
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Fig. 4-2 - Typical Line Printer Output 
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4.6 SUMMARY 

The ma.in handicap in using this program is the long run time required 
(approximately 20 minutes, for 0.5 second of real time). Therefore, to simu- 
late a 20-second run would take over 13 hours of 1108 central processor time. 
However, short runs have shown some interesting phenomena, such as the 
large variations that the random particle distributions make on the signal 
amplitude in the various filters. More runs using the program are desirable. 
Possibly some program changes can be made to speed up the program. 
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Appendix 

ACCEPTANCE TEST PROCEDURE 


A.1 ELEVATION 
A, 1,1 Automatic Mode 
A.i.1,1 Upper Scan Limit Set 

For test purposes the scan limit is defined as being the angle deter- 
mined by the voltage output of the D/A converter for the command word when 
the latter voltage is at its maximum positive value. The maximum angle is 
given by: 


9 


max 


10.24 (degrees) 


(A.1) 


The value of a-s determined by (A.1) shall be functionally related 

to the angle set on the upper thumbwheel in accordance with Table A-1. 
Table A-1 gives the effect of the 5 bit truncation used in this comparator, 

A, 1.1.2 Automatic Mode — Lower Scan Limit Set 

For test purposes the scan limit is defined as being the angle deter- 
mined by the voltage output of the command word D/A converter when the 
latter is at its minimum value. The minimum angle is given by 


9 = 10,24 V . (degrees) 

max min ' ® 


(A, 2) 


The value at is given by (A. 2) and shall be functionally related 

to the lower limit thumbwheel setting as given in Table A- 2. 
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Table A-1 

ELEVATION SCAN UPPER LIMIT 


Actual Scan 
Upper Limit 
(Digital) 

Angle 

Value 

Thumbwheel Setting 
(Upper Limit) 


10010 

57.6 

59, 58 

Note 1 

10001 

54.4 

57, 56, 55 


10000 

51.2 

54, 53, 52 


01111 

48.0 

51, 50, 49, 48 


OHIO 

44.8 

47, 46, 45 


01101 

41J> 

44, 43, 42 


01100 

38,4 

41, 40, 39 


01011 

35.2 

00 


01010 

32.0 

35, 34. 33, 32 


01001 

28.8 

31, 30. 29 


01000 

25,6 

28, 27, 26 


00111 

22.4 

25, 24, 23 


00110 

19.2 

22, 21, 20 


00101 

16.0 

19, 18, 17, 16 


00100 

12.8 

15, 14, 13 


00011 

9.6 

12, 11, 10 

Note 2 


Note i: Mechanical stops limit thumbwheel to 59 deg maximum 
Note Z: Mechanical stops limit thumbwheel to 10 deg minimum 
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Table A-2 

ELEVATION SCAN LOWER LIMIT 


Actual Scan 
Lower Limit 
(Digital) 

— 

0 

max 

Thumbwheel Setting 
(Lower Limit) 


01111 

51.175 

49,48 

Note 1 

OHIO 

47.97 5 

47, 46, 45 


01101 

44.775 

44, 43. 42 


01100 

41.575 

41,40, 39 


01011 

38.375 

38, 37. 36 


01010 

35.175 

35, 34, 33, 32 


01001 

31.97 5 

31, 30, 29 


01000 

28.775 

28, 27, 26 


00111 

25.57 5 

25, 24, 23 


00110 ' 

22.375 

22, 21, 20 


00101 

19.175 

19, 18, 17, 16 


00100 

15.975 

15, 14, 13 


00011 

12.775 

12, 11, 10 


00010 

9.57 5 

9, 8, 7 


00001 

6.375 

6. 5, 4 


00000 

3.175 

3 or smaller 



Note 1: Table applies so long as the difference between upper 

and lower limits on the thumbwheel is at least 10 deg; 
if not, then - 9*6 which is accomplished 

automatically. 
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A. 1.1, 3 Limit Tests and Procedure — Automatic Mode 
Test E-1 


The output at the command word D/A converter shall be connected to 
a calibrated strip chart recorder with a convenient scale factor. The control 
settings should be as follows: 

a. Mode — automatic 

b. Frequency — 0.1 Hz 

c . Lower limit — 0 deg 

Starting at 10 deg, the upper limit shall be increased manually by 1 deg 
at the completion of at least one complete scan cycle. The values of 0 as 
read from the recorded data {using Eq. (A.l)) shall conform to Table A-1 within 
+2 deg. 

Test E-2 

Using the same test set-up, the controls should be set as follows; 

a. Mode — automatic, normal 

b. Frequency — 0.1 Hz 

c. Upper limit 59 deg 

Starting at 0 deg, the lower limit should be increased manually by 

1 deg at the completion of at least one scan cycle. The values of 0 . 

mm 

should conform to Table A-2 within +2 deg. {Note; the difference between 
upper and lower limits should not be made less than 10 deg.) 
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Test E-3 — Scan Frequency 


For each of the conditions given in Table A-3 below, manually program 
the scan frequency through the five available values, and determine the scan 
frequency from the strip chart recording, using time calibration marks or 
known chart speed. The frequency should correspond to the dial setting to 
within +10 percent. 

Table A-3 


Test 

Condition 

Upper 

Limit 

(deg) 

Lower 

Limit 

(deg) 

Span 

(deg) 

Center 

(deg) 

1 

59 

44 

15 

51.5 

2 

39 

24 

15 

31.5 

3 

19 

4 

15 

11.5 

4 

59 

29 

30 

44.0 

5 

34 

4 

30 

19.0 

6 

59 

9 

30 

20.0 


Test E-4 — Crossover 


a. Control settings: automatic, manual 

b. Frequency: 0.1 Hz 

With the same test set-up, set the upper limit successively to 20, 40, 

and 59 deg. At each setting, move the lower limit toward the upper limit; 

the lower limit on the thumbwheel should be automatically replaced by =9.6 

U L 

deg whenever the difference is less than 10 deg. 

A. 1.2 Manual Mode 

The command voltage should be recorder monitored as in the previous 

tests. 
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Test E-5 

Depression of the up/down slewing switches should produce monotone 
variation of the command voltage between the limits of 89.6 and 3,2 deg in the 
selected direction. The indicated read-out angle should correspond to the 
chart recording to within +0.20 deg. The slew speeds should be 34 deg/sec 
and 3,4 deg/sec to within J^10%, 

Depression of the single increment button, when the speed setting is 
"fast," should enable single increments of no greater than 0.2 deg in either 
selected direction. 


A. 1.3 Sector Scan Mode 

Tests of this mode shall employ the same recorder test set up as in 
the previous tests 


a. Mode — automatic, sector 

b. Frequency — 0.1 Hz 


Test E-6 


Press the "clear” button to index the control counters. 


a. Depress the "span increase" lever for several seconds to move 
the counter control word to maximum. 

b. Depress the "span diminish" lever for several seconds to run 
the span counter down to its automatic limit of 3.175 deg. 

c. Depress the "center location" lever in the up and down directions; 
the upper and lower limits as observed on the recorder should 
move in unison (within ^3.2 deg) maintaining a span of 6.4 deg. 

d. With the center of the scan excursion stopped, depression of the 
span control in the selected direction should provide monotone 
expansion or contraction of the scan excursion about the original 
center angle. 
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e. The span should be capable of expansion to +30 deg maximum and 
the center should be capable of excursion from 1.0 to 50 deg. The 
upper limit should be automatically restrained from falling below 
3.2 deg (these will inhibit symmetry of the scan limits about the 
center line under the appropriate circumstances). 


A. 2 RANGE 

The definitions and tests for range are identical to those for elevation, 
with the following exceptions. 

a. Range values correspond to ten times the counterparts in angle, 

b. Tests specifying 0.1 Hz scan frequency for 0 should be run at 
1 Hz for range. Due to the higher speeds, changes in limit 
settings should be made every 5 seconds rather than once per 
cycle. 

c. The frequency tests of test E-3 should be run at 0.1, 1, 2, 3,4 
and 5 Hz. 


A. 3 SERVO PERFORMANCE 
A. 3.1 Elevation Servo 
A. 3. 1.1 Static Accuracy Manual Mode 
Test E-7 


With the elevation system in the manual mode, the angle shall be slewed 
manually in 5-deg increments from 5 to 55 deg, as indicated on the LED read- 
outs, and returned down through the same values. The angle indicated on the 
vernier dial should correspond to the LED reading to +0.2%. 


A. 3. 1.2 Dynamic Response — Automatic Mode 

At the highest scan rate (0.5 Hz) and maximum excursion (*^55 deg) the 
phase shift between driving sawtooth command and mirror motion should not 
exceed 20 deg. The scan excursion should be within 5 deg of the command 
angle at both turn-around points. The period during which the scan velocity 
falls in absolute value to less than 75% of the command velocity shall not 
occupy more than 10% of the scan period. 
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These values shall be determined by the strip chart recording of the 
command voltage and feedback pot voltage after the latter has been aligned 
and the gains and offsets adjusted. 

A. 3. 2 Range Servo 

A. 3. 2.1 Static Accuracy 

Since the actual location of the focus (point of maximum Intensity of the 
illumination) is not observable within the control system, static accuracy 
tests should encompass at least one determination of the correspondence be- 
tween the secondary mirror position and the location of the actual physical 
focus, using the following technique. 


Test R-8 


a. A Doppler wheel target shall be setup externally at range R^; 
the range servo shall be manually set to give an LED readout 
corresponding to R^, and, with the loop closed and elevation 
angle correctly aimed, the secondary mirror platform position 
is adjusted to give maximum voltage from a fixed I.F. bandwidth 
(non -scanning) spectrum analyzer signal power output terminal. 

b. The target is moved to range R 2 and the spectrum analyzer output 
is again maximized by manually slewing the range to whatever 
new value is required to accomplish this. Both the LED readout 
and the optical encoder D/A converter voltage output correspond- 
ing to this maximum should be recorded. The latter value (after 
multiplication by the scale factor relating encoder angle to mirror 
translation) is to be taken as the mirror position that the positions 
should satisfy. 


(di 


>z> = ^ - tj 


where 

dj = mirror position 
F = system focal length 
R = target range 


For this test, Rj^ = 100 meters, ^2 ~ meters. 
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Test R-9 

The output of the optical encoder D/A converter (which measures the 
secondary mirror position) shall be measured with a precision digital volt- 
meter while the range is slewed to successive positions 50 meters apart, 
from 50 meters to 550 meters, and back to 50 meters stopping at the same 
values and maintaining monotone directional sense. The range values are 
those indicated by the LED range readouts. After conversion by the scale 
factor, the mirror position values should satisfy. 

d. = F^/R. 

(within an additive constant) to within 1% of the corresponding values based 
on the Rj^ observed on the LED, At each position of static range, the re- 
constructed range voltage shall also be measured, and after scale factor 
conversion, this should be within 1% of the command range, 

A. 3. 2. 2 Dynamic Performance 

Test R-10 

The command range input voltage and D/A (mirror position) voltage 
shall be simultaneously recorded on strip chart. The excursion shall be set 
to span 6.4 to 480 meters. A run shall be made at 0.1 Hz scan frequency 
for comparison purposes. When the scan frequency is increased to 5 Hz, 
the phase shift (between input sawtooth minima and mirror position maxima) 
shall not exceed 10 deg, and the excursion shall not depart from the low fre- 
quencies by more than 5%. 



